arXiv:1504.07050vl [gr-qc] 27 Apr 2015 


Cosmology and Gravitation: the grand scheme for High-Energy Physics 


P. Binetruy 

Universite Paris Diderot, Paris, France 

Abstract 

These lectures describe how the Standard Model of cosmology (ACDM) has 
developped, based on observational facts but also on ideas formed in the con¬ 
text of the theory of fundamental interactions, both gravitational and non- 
gravitational, the latter being described by the Standard Model of high en¬ 
ergy physics. It focuses on the latest developments, in particular the precise 
knowledge of the early Universe provided by the observation of the Cosmic 
Microwave Background and the discovery of the present acceleration of the 
expansion of the Universe. While insisting on the successes of the Standard 
Model of cosmology, we will stress that it rests on three pillars which involve 
many open questions: the theory of inflation, the nature of dark matter and 
of dark energy. We will devote one chapter to each of these issues, describ¬ 
ing in particular how this impacts our views on the theory of fundamental 
interactions. More technical parts are given in italics. They may be skipped 
altogether. 

1 A not so brief history of modern cosmology 

Cosmology has been an enquiry of the human kind probably since the dawn of humanity. Modern 
cosmology was born in the early XX th century with the bold move of Einstein and contemporaries to 
apply the equations of general relativity, the theory of gravity, to the whole Universe. This has led to many 
successes and/or suprises, the most notable of which being presumably the discovery of extra-galactic 
objects which recede from our own Galaxy, i.e. the discovery of the expansion of the Universe 00 - 
This led to the development of the Big Bang theory, with the early Universe being a hot and dense 
medium (a prediction confirmed by the discovery of the cosmic microwave background by Penzias and 
Wilson in 1965 [3j), and thus a laboratory for studying elementary particles. A picture thus emerged in 
the 1970s, not only based on the theory of gravity, but also on non-gravitational interactions described by 
the Standard Model of high energy physics, which was being finalized at the same time (its experimental 
confirmation would take another 40 years and has culminated in the discovery of the Fliggs particle in 
2012 ). 

A first success of the particle physics approach to cosmology has been the understanding of the 
abundancy of light elements in the 80s. This was the first quantitative success of cosmology. Mean¬ 
while, the development of gauge symmetries and the understanding of the role of spontaneous symmetry 
breaking in fundamental interactions led the community to focus its attention on phase transitions in the 
early Universe, in particular associated with the quark-gluon transition, the breaking of the electroweak 
symmetry or even of the grand unified symmetry. It is in this context that, in the early 80s, the theory of 
inflation was proposed @0 to solve some of the mysteries of the standard Big Bang theory. 

The theory of inflation included a model for the genesis of density fluctuations responsible for 
the formation of large scale structures, such as galaxies or clusters of galaxies: the quantum fluctuations 
during the exponential (de Sitter) expansion. But this implied the presence of fluctuations in the otherwise 
homogeneous and istropic Cosmic Microwave Background. Such fluctuations were observed by the 
COBE satellite, at the level of one part in 100 000. Generic models of inflation predicted also in a very 
elegant manner that space (not spacetime!) is flat, any spatial curvature being erased by the exponential 
expansion. According to Einstein’s equations, this implied that the average energy density in the Universe 
had the critical value p c r\j 10 26 kg/m 3 . 


This was a prediction not supported by observation. It was known since the 1930s that there was a 
significant amont of non-luminous -or dark- matter in the Universe: in 1933, Fritz Zwicky, by studying 
the velocity distribution of galaxies in the Coma cluster, had identified that there was 400 times more 
mass than expected from their luminosity. This had been confirmed by studying subsequently the rotation 
curves of many other galaxies. But the total of luminous and dark matter could not account for more than 
30% of the critical energy density (other components like radiation are subdominant at present times). 
Models of open inflation were even constructed to reconcile inflation with observation. 

The clue came in 1999 |6j[7j when it was observed that the expansion of the Universe is presently 
accelerating. Since matter or radiation tend to decelerate the expansion, one has to resort to a new form 
of energy, named dark energy, to understand this acceleration. Was this the component which would 
provide the missing 70% to account for a total energy density p c and thus a spatially flat Universe? 
The answer came from a more precise study of the fluctuations in the CMB through the space mission 
WMAP (and more recently Planck): they conclude indeed that these fluctuations are consistent with 
spatial flatness. 

The latest cosmology results from the Planck mission, released this year, have confirmed the 
predictions of the simplest models of inflation, a rather remarkable feat since they are associated with 
dynamics active in the first fractions of seconds after the big bang, and they allow tu fully understand the 
imprints observed 350 000 years after the big bang. 

We thus have at our disposal a Standard Model of cosmology which is sometimes compared with 
the Standard Model of high energy physics: in both cases, no major experimental/observational data 
seems to be in conflict with the Model. There is however one big difference. The Standard Model of 
cosmology rests on three “pillars” -inflation, dark matter, dark energy- which are very poorly known: 
we have at present no microscopic theory of inflation, and we ignore the exact nature of dark matter or 
dark energy. 

For example, there are convincing arguments that dark matter is made of weakly interacting mas¬ 
sive particles of a new type, and a large experimental programme has been set up to identify them. Their 
discovery would be of utmost importance because this would be the first sign of physics beyond the 
Standard Model of particle physics. But it remains a possibility -though not a favored one- to explain 
the observed facts through a modification of gravity at different scales (from galaxies to clusters and 
cosmological scales). Finally, axion dark matter would be a minimal extension fof the Standard Model, 
accounting for dark matter. 

The discovery of the Higgs has provided us with the first example of a fundamental scalar field 
(at least fundamental at the scale where we observe it), this a welcome for cosmology since micro¬ 
scopic models of accleration of the expansion of the Universe -whether inflation or dark energy- make 
heavy use of such fields. They have the double advantage of being non-vanishing without breaking the 
symmetries of spacetime (like Lorentz symmetry) and of having the potential of providing an unclus¬ 
tered background, they also appear naturally in the context of extensions of the Standard Model, like 
supersymmetry or extra dimensions. 

Scalar fields thus provide valuable toy models of inflation or dark energy. However, such toys 
models are difficult to implement into realistic high energy physics models because of the constraints 
existing on physics beyond the Standard Models of particle physics and cosmology. 

In what follows, we will review the successes of the Standard Model of cosmology, focussing on 
the most recent results. And we will consider the three pillars of this Model —inflation, dark matter, 
dark energy- and identify the most pressing questions concerning these three concepts. 

In what follows some more technical discussions are given in italics. They can be skipped alto¬ 
gether. Background material is also provided in Appendices, most of it being intended for the reader who 
wants to reproduce some of the more advanced results. 

We start in this first Chapter with an introduction, partly historical, to the main concepts of cos- 


2 


mology. 


1.1 Gravity governs the evolution of the Universe 

The evolution of the universe at large is governed by gravity, and thus described by Einstein’s equations. 
We recall that general relativity is based on the assumption made by Einstein that observations made 
in an accelerating reference are indistinguishable from those made in a gravitational field (as illustrated 
on a simple example in Fig. |T]). This has several consequences, the most notable of which being that 


constant velocity 


acceleration lg 



Fig. 1: (a) When the astronaut in a rocket in uniform motion (constant velocity) drops the balls, they float around, 
(b) If the rocket accelerates by 1 g, the balls fall to the ground and the astronaut has no way to identify whether 
this due to acceleration or to a gravitational field fi.e. the Earth attraction). 


the curves followed by light (null geodesics) are not straight lines (see Fig. [2]). Einstein’s equations are 
highly non-linear second order differential equations for the metric g^^\ They reacQ 

Gmi = R^,v 2 Sfj.i'R — j^T^ii T A g^u • (3) 


where I! IU , is the Ricci tensor, R the associated curvature scalar (see Appendix [fijl and T /u/ the energy- 
momentum tensor; finally A is the cosmological constant which has the dimension of an inverse length 
squared. 

1 We recall that the invariant infinitesimal spacetime interval reads 

ds 2 = g^dx^dx 1 ' . (1) 


The metric signature we adopt throughout is Einstein’s choice: (+, —, —, —). 

~ These equations can be obtained from the Einstein-Hilbert action: 

s -s 

where the generic fields if contribute to the energy-momentum: T fv = (2/g)(6Sm/Sg IJ ' 1 ') 


--H-X 


+ Sm (t/ 1 , giin) , 


( 2 ) 
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constant velocity 


acceleration 



Fig. 2: If light propagates along straight lines in the rocket in uniform motion (a), it must deviate from a straight 
line when the rocket accelerates (b). Thus it must propagate along curved lines in gravitational fields. 


Thus Einstein’s equations relate the geometry of spacetime (the left-hand side of ([3])) with its 
matter field content (the right-hand side). 

Einstein equations are field equations. Of which field? This is better understood in the weak gravi¬ 
tational field limit, that is in the limit of an almost flat spacetime. In this case, the metric is approximated 
by 

9nv ~ Vuv + hpv(x) , (4) 

where r) jW is the Minkowski metric and h tw {x) is interpreted as the spin-2 graviton field. In particular, 
we note that 

9oo = 1 + 2<3? , (5) 

where <b is the Newtonian potential which satisfies the Poisson equation A4> = AttG N p. 

Einstein used his equations not just to describe a given gravitational system like a planet or a star 
but the evolution of the whole universe. In his days, this was a bold move: it should be remembered 
how little of the universe was known at the time these equations were written. “In 1917, the world was 
supposed to consist of our galaxy and presumably a void beyond. The Andromeda nebula had not yet 
been certified to lie beyond the Milky Way.”[Pais |J8| p. 286] Indeed, it is in this context that Einstein 
introduced the cosmological constant in order to have a static solution (until it was observed by Hubble 
that the universe is expanding) for the universe. 

More precisely ||8j, Einstein first noticed that a slight modification of the Poisson equation, namely 

A<J> — A4> = 4-7T G N p (6) 

allowed a solution with a constant density p (4> = —AttG n p/X) and thus a static Newtonian universe. In 
the context of general relativity, he found a static solution of ([3]) under the condition that 

A = -^ = 4vr G N p , (7) 
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where r is the spatial curvature (see Exercise 1-1). It was soon shown that this Einstein universe is 
unstable to small perturbations. 

Exercise 1-1 : Consider the following metric 


5oo = l, 9ij = -Sij + J'* 3 r2 , x 2 = J2 x i ■ W 

2—1 

a) Show that it is a solution of Einstein’s equations (jd]j in the case of non-relativistic matter with 
a constant energy density p satisfying the condition 0 

b) Prove that, in the Newtonian limit, one recovers 0 

Hints: a) T l jk = r -2 [xt^jk ~ XiXjXk/(x 2 — r 2 )] which gives Rij = —2 gij/r 2 . 
b) hi the Newtonian limit, Goo Apoo wdh poo given by 0 


1.2 An expanding Universe 

Other solutions to the Einstein equations were soon discovered. The first exact non-trivial one was 
found in late 1915 by Schwarzschild, who was then fighting in the German army, within a month of the 
publication of Einstein’s theory and presented on his behalf by Einstein at the Prussian Academy in the 
first days of 1916 [9j, just before Schwarzschild death from a illness contracted at the front. It describes 
static isotropic regions of empty spacetime (A = 0), such as the ones encountered in the exterior of a 
static star of mass M and radius R: 

dg2 = ^ _ 2G^ M ^ df2 _ ^ _ 2G^M\ _1 dr 2 _ r 2 d0 2 _ r 2 sin 2 ed(j) 2 (9) 

In 1917, de Sitter proposed a time-dependent vacuum solution in the case where A / 0 E0HD: 

ds 2 = dt 2 - e 2Ht dy? , H 2 = A 2 /3 . (10) 


But, since there exists time-dependent solutions, why should the Universe be static? This led to the so- 
called “Great Debate” between Harlow Shapley and Heber D. Curtis in 192(0 H. Shapley was supporting 
the view that the Universe is composed of only one big Galaxy, the spiral nebulae being just nearby gas 
clouds (he was also arguing -rightfully- that our Sun is far from the center of this big Galaxy). H.D. 
Curtis, on the other hand considered, that the Universe is made of many galaxies like our own, and 
that some of these galaxies had already been identified, in the form of the spiral nebulae (and he was 
supporting the view that our Sun is close to the centre of our relatively small Galaxy). 


In 1925, Edwin Hubble studies, with the 100 inch Hooker telescope of Mount Wilson, the Cepheids 
which are variable stars in the Andromeda nebula M31. He shows that the distance is even greater than 
the size proposed by Shapley for our Milky Way: M31 is a galaxy of its own, the Andromeda galaxy (at 
a distance of 3.10 6 light-years from us) 121. 

In 1929, Hubble (2) discovers, by combining spectroscopic measurements with measures of dis¬ 
tance, that galaxies at a distance d from us recede at a velocity v following the law: 


v ~ Hod , (11) 

the constant of proportionality being henceforth called the Hubble constant (see Fig. [30 As a conse¬ 
quence of the Hubble law, the Universe is expanding! 

3 see http://apod.nasa.gov/diamondjubilee/debate20.html 

4 Hubble’s result was actually anticipated by G. Lemaitre who published the law ® two > ears earlier 0- 
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Fig. 3: Hubble plot (radial velocity vs distance) of extragalactic “nebulae” obtained by Hubble in his original 
article of 1929 |2| and by Hubble and Humason two years later fl3] . 


The velocities v are measured by Hubble through the Doppler shift of spectroscopic lines of the 
light emitted by the galaxy: if A em it is the wavelength of some spectrocopic line of the light emitted by 
a galaxy (receding from us at velocity v) and A 0 b s of the corresponding line in the light observed by us, 
then 


l + z = 


Aobs 


*emit 


1 + v/c 

y/l — V 2 /c 2 


~ 1 + v/c . 


( 12 ) 


Hence v ~ zc. 


Let us note that stars within galaxies, such are our Sun, are not individually subject to the expan¬ 
sion: they have fallen into the gravitational potential ofthe galaxy and are thus not receding from one 
another (see Exercise 1-2). This is why one had first to identify extragalactic objects before discovering 
the expansion: only galaxies, or clusters of galaxies, recede from one another. 


Exercise 1-2 : The purpose of this exercise is to show that, in the case where a fluctuation of 
density appears (such as when galaxies form), the massive objects (stars) decouple from the general 
expansion to fall into the local gravitational potential /|/4||/5[ /. 

In a matter-dominated universe of uniform density p, a perturbation appears in the form of a 
sphere with uniform excess density A p. Assuming that the gravitational field inside the sphere is de¬ 
scribed by the Robertson-Walker metric with positive curvature constant A k > 0, the evolution is gov¬ 
erned by the equations 

a 2 A k 1. . 

H- yy = x(P + Ap)) (13) 

a z cr 6 

a 3 (p + Ap) = C. (14) 

where C is a constant and we neglect the cosmological constant. 
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a) Show that the solution is given parametrically by 


a{t) 


t 


C 

6Ak 

to + 


(1 — cos rj) , 
C 


6 A/c 3 / 2 


(77 — sin r/) 


(15) 

(16) 


Note that this does not assume that A p is small. 

b) Show that, as t —»• to, one has a(t) = (C/ 3 ) 1//3 [3 (t — to)/2 ] 2 ^ 3 ~ (t — to ) 2 / 3 , as in the rest of 
the matter-dominated universe. Verify that, whereas p = (4/3)(t —to) 2 (~ a~ 3 ), A p = (9/5)A ka~ 2 ~ 
(t - to)- 4/3 . 

c) t/mv Zong does it take before the system starts to collapse to a bound system or to a singularity? 


Hints: a) ad 2 + aAk = C/3 

b) Make an expansion in i] to the subleading order (q is small when t 


to)- 


c) From (15) the expansion stops and the collapse starts at ij = tt or t — to = 7rC'/(6AA:' !/ 2 ). 


Exercise 1-3 : Identify the redefinition of coordinates t = t(t, r), r = r(f, r) which allows to write 
the de Sitter metric (jTO]) 

(IV) 


ds z = dt z - e 2Ht {dr 2 + rW) , 


into the following form: 


^2 


R 2 , 


ds z = I 1 — — 7j- ) dt 2 — ( 1 — —y ) dr z — r z dVt z , Rh = H 


^2 \ -t 


R 2 h 


j^2 -2 jo 2 


- 1 


( 18 ) 


The first equation is known as the flat form of the de Sitter metric (see ( |T9| ) below), the second one is the 
static form (compare with the Schwarzschild metric Q). 


Hints: e Ht = e Ht y/l — H 2 r 2 and r = re Ht . 


1.3 Friedmann-Lemaitre-Robertson-Walker universe 

As one gets to larger and larger distances, the Universe becomes more homogeneous and isotropic. Under 
the assumption that it reaches homogeneity and isotropy on scales of order 100 Mpc (1 pc = 3.262 
light-year = 3.086 X 10 16 m) and larger, one may first try to find a homogeneous and isotropic metric 
as a solution of Einstein’s equations. The most general ansatz is, up to coordinate redefinitions, the 
Robertson-Walker metric: 


ds 2 
'yijdx l dx J 


c 2 dt 2 — a 2 (t ) 7 ijdx l dxf 
1 ^ kr 2 + r2 (^ 2 + Sin2 0 # 2 ) ) 


(19) 

( 20 ) 


where a(t) is the cosmic scale factor, which is time-dependent in an expanding or contracting universe. 
Such a universe is called a Friedmann-Lemaitre universe. The constant k which appears in the spatial 
metric 7 ^ can take the values ±1 or 0: the value 0 corresponds to flat space, i.e. usual Minkowski 
spacetime; the value +1 to closed space (r 2 < 1) and the value —1 to open space. Note that r is 
dimensionless whereas a has the dimension of a length. From now on, we set c = 1, except otherwise 
stated. 

For the energy-momentum tensor that appears on the right-hand side of Einstein’s equations, we 
follow our assumption of homogeneity and isotropy and assimilate the content of the Universe to a perfect 
fluid: 

T^ v = -pg^v + (p + P)U,jU v , (21) 
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where U 1 ' is the velocity 4-vector (JJ t = 1, U l = 0). It follows from ( |2lj ) that T tt = p and T rJ = a' 2 pj r] . 
The pressure p and energy density p usually satisfy the equation of state: 


p = wp 


( 22 ) 


The constant w, called the equation of state parameter, takes the value w ~ 0 for non-relativistic matter 
(negligible pressure) and w = 1/3 for relativistic matter (radiation). In all generality, the perfect fluid 
consists of several components with different values of w. 

One obtains from the (0, 0) and (i, j) components of the Einstein equations <[3]> (see Exercise B-l 
of Appendix [77]) : 


cr k 

— H—y ) = 87 tG n p + A, 

a z a z 


a 2 + 2ad + k = — 8irG N a A p + a z A, 


(23) 

(24) 


where we use standard notations: a is the first time derivative of the cosmic scale factor, a the second 
time derivative. 

The first of the preceding equations can be written as the Friedmann equation, which gives an 
expression for the Hubble parameter H = a/a measuring the rate of the expansion of the Universe: 


H 2 — yy — - (A + 87 tG n p) - . 


a z 


a z 


(25) 


The cosmological constant appears as a constant contribution to the Hubble parameter. We note that, 
setting k = 0 and p = 0, one recovers de Sitter solution (10): a 2 /a 2 = A/3 i.e. aft ) ~ e Ht . For the time 


being, we will set A to zero and return to it in subsequent chapters. 
Next, we note that, assuming k = 0, we have at present time: 


3 H 2 _ 

P = 8ttGZ = Pc 


10 _ 2 b kg/nU 


(26) 


where Hq is the Hubble constant, i.e. the present value of the Hubble parameter. This corresponds to 
approximately one galaxy per Mpc 3 or 5 protons per m 3 . In fundamental units where h = c = 1, this 
of the order of (lO~ 3 eV) 4 . We easily deduce from (25 (that space is open (resp. closed) if at present 


is 


time p < p c (resp. p < p c ). Hence the name critical density for p c . 


Friedmann equation (25 l can be understood on very simple grounds: since the universe at large 
scale is homogeneous and isotropic, there is no specific location and motion in the universe should not 
allow to identify any such location. This implies that the most general motion has the form v(f) = H (t)x 
where x and v denote the position and the velocity and H ( t) is an arbitrary function of time. Since 
v = x, one obtains x = a(f)r, where r is a constant for a given body (called the comoving coordinate) 
and aft) is related to Hft) through H = a/a. Now, consider a particle of mass m located at position x: 
the sum of its kinetic and gravitational potential energy is constant. Denoting by p the energy density of 
the (homogeneous) universe, we have 


1 o 47T . ,2 

-mv-— Ec N mp |x| = cst. 

Z O 


Writing this constant —kmr z /2, we obtain from (27) 


a 2 8n k 

a? = T &Np ~a 2 


(27) 


(28) 


which is nothing but Friedmann equation (251 (with vanishing cosmological constant). 








Friedmann equation should be supplemented by the conservation of the energy-momentum tensor 
which simply yields: 

P= -3H{p + p) . (29) 


Hence a component with equation of state (22 1 has its energy density scaling as p ~ a (t)“3(i+“). Thus 
non-relativistic matter (often referred to as matter) energy density scales as a In other words, the 
energy density of matter evolves in such a way that pa 3 remains constant. Radiation scales as a 1 and a 
component with equation of state p = — p (w = — 1) has constant energy density^] 

We note for future use that, if a component with equation of state ( [22] ) dominates the energy density 
of the universe (as well as the curvature term —k/a 2 ), then (281 has a scaling solution 

2 


a(t) ~ t v 


with v = 


(30) 


3(1 + w) 

For example, in a matter-dominated universe, a(t) ~ f 2 / 3 , and in a radiation-dominated universe, a(t) ~ 

t 1 / 2 . 

Differentiating the Friedmann equation with respect to time, and using the energy-momentum 
conservation (|29|), one easily obtains 


a = -—a(3p + p) + a- . 


(31) 


This allows to recover (24 1 from Friedmann equation and energy-momentum conservation. 


1.4 Redshift 

In an expanding or contracting universe, the Doppler frequency shift undergone by the light emitted from 
a distant source gives a direct information on the time dependence of the cosmic scale factor a(t). To 

obtain the explicit relation, we consider a photon propagating in a fixed direction {9 and <fi fixed). Its 

equation of motion is given as in special relativity by setting ds 2 = 0 in ( fl9] ): 

c 2 dt 2 = a 2 (t) dr . (32) 

1 — kr z 

Thus, if a photon (an electromagnetic wave) leaves at time t a galaxy located at distance r from us, it 
will reach us at time to such that 

f = r dr 03) 

Jt a (t) Jo Vl — hr 2 

The electromagnetic wave is emitted with the same amplitude at a time t + T where the period T is 
related to the wavelength of the emitted wave A by the relation A = cT. It is thus received with the same 
amplitude at the time to + To given by 

rto+To cdt _ [ r dr 

it+T a(t) Jo ' 

the wavelength of the received wave being simply Ao = cTo- Since To,T <C to,t, we obtain from 
comparing 

rTn rT 1 An ao 

(35) 


f 

Jt 4 


cTo cT 
ao a(t ) 


i.e. 


Ap _ _ 

A a(t) 


where ao is the present value of the cosmic scale factor. 

Defining the redshift parameter z as the fractional increase in wavelength z = (Ao — A)/A, we 

have 


1 + 2 = 


a o 
a(t) 


(36) 


One may thus replace time by redshift since time decreases monotonically as redshift increases. 


5 The latter case corresponds to a cosmological constant as can be seen from 1 23[24[ where the cosmological constant can 
be replaced by a component with pA = — Pa = A/(87rG N ). 
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1.5 The universe today: energy budget 

The Friedmann equation 

H 2 = \ (A + 8v tG nP ) - 4 • (37) 

or 3 or 

allows to define the Hubble constant Hq, i.e. the present value of the Hubble parameter, which sets the 
scale of our Universe at present time. Because of the troubled history of the measurement of the Hubble 
constant, it has become customary to express it in units of 100 km.s _1 .Mpc 1 which gives its order of 
magnitude. Present measurements give 


ho 


II 

100 km.s -1 .Mpc _1 


0.7 ±0.1 


The corresponding length and time scales are: 

l Ho = 4r = 3000 hg - 1 Mpc = 9.25 x 10 25 h^ 1 m, 
-no 

t H o = — = 3.1 x 10 17 JiQ 1 s = 9.8 Hq 1 Gyr. 
no 


(38) 

(39) 


It has become customary to normalize the different forms of energy density in the present Uni¬ 
verse in terms of the critical density p c = 3Hq/( 8irG N ) defined in (261. Separating the energy density 
p M0 presently stored in non-relativistic matter (baryons, neutrinos, dark matter,...) from the density p R0 
presently stored in radiation (photons, relativistic neutrino if any), one defines: 


^M = 


n R = ^, n A = 

Pc Pc 


A 


3 H 2 ' 


Vk — 


k 


a l H l ' 


(40) 


The last term comes from the spatial curvature and is not strictly speaking a contribution to the energy 
density. One may add other components: we will refrain to do so in this Chapter and defer this to the last 
one. 

Then the Friedmann equation taken at time to simply reads 


(7 m + (7 fl + AIa — 1 — Ufc . 


(41) 


Since matter dominates over radiation in the present Universe, we may neglect Q lt in the preceding 
equation. As we will see in the next Chapters, present observational data tend to favor the following set 
of values: Ufc ~ 0 (see Section |2T| and Cl M ~ 0.3, ~ 0.7 (see Section [5T] and Fig. 22 1 . The matter 

density is not consistent with the observed density of luminous matter (Hi uminous ~ 0.04) and thus calls 
for a nonluminous form of matter, dark matter which will be studied in Section |4| 

Using the dependence of the different components with the scale factor a(t) = «q/(1 + z), one 
may then rewrite the Friedmann equation at any time as: 


H 2 (t) 
or H 2 (z) 


m 


ft a + U. 


_ao_ 

a(t ) 


+ ru 


Qq 

a(t) 


+ S7i 


Qq 

a(t) 


H 2 [(7 m (1 + ^) 3 + (7 fl (l + z ) 4 + (7fc(l + z) 2 + (7a] ■ 


(42) 

(43) 


where ao is the present value of the cosmic scale factor and all time dependences (or alternatively redshift 
dependence) have been written explicitly. We note that, even if Q lt is negligible in ( [4T] ), this is not so 
in the early Universe because the radiation term increases faster than the matter term in ( |42| ) as one gets 
back in time (i.e. as a(t) decreases). If we add an extra component X with equation of state p x = w x p x , 
it contributes an extra term Q x (ao/a(t)) 3 ^ 1+Wx ' > where U x = p x /p c - 
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An important information about the evolution of the universe at a given time is whether its ex¬ 
pansion is accelerating or decelerating. The acceleration of our universe is usually measured by the 
deceleration parameter q which is defined as: 


da 



(44) 


Using ([37]) of Section[2]and separating again matter and radiation, we may write it at present time to as: 


qo = 


R 2 

n o 


t=t 0 


= 

2 3 


+ Qr — 


(45) 


Once again, the radiation term Q ;( can be neglected in this relation. We see that in order to have an 
acceleration of the expansion (qo < 0), we need the cosmological constant to dominate over the other 
terms. 


We can also write the deceleration parameter in (44 1 in terms of redshift as in (43) 


H o 

2H(z) 2 


[( 2 M (1 + 2 ) 3 + + z) 4 — 2f2 A ] 


(46) 


This shows that the universe starts accelerating at redshift values 1 + z (2U a Mw) 1/3 (neglecting 

£l R ), that is typically redshifts of order 1. 

The measurements of the Hubble constant and of the deceleration parameter at present time allow 
to obtain the behaviour of the cosmic scale factor in the last stages of the evolution of the universe: 


a(t) = a 0 


1 + 


t ~ t 0 

t H 0 


qo (t -1 0 ) 2 


2 


t 


2 

Ho 


+ ■■■ 


(47) 


1.6 The early universe 

Table [T] summarizes the history of the universe in the context of the big bang model with an inflationary 
epoch. We are referring to the different stages using time, redshift or temperature. The last two can be 
related using the conservation of entropy. 

Indeed, one can show, using the second law of thermodynamics (TdS = dE + pdV) that the 
entropy per unit volume is simply the quantity 


_ S _ p + p 
S ~ V T 


(48) 


and that the entropy in a covolume sa 3 remains constant. The entropy density is dominated by relativistic 
particles and reads 


» = \ds{T)a BB T 3 , (49) 

9 s ( T ) = £ + ^ £ 9, (|) 3 , (50) 

bosons i x / fermions i x / 

where the sum extends only to the species in thermal equilibrium and a BB = tt 2 /,' 4 /(loc 3 /; 3 ) = 7.56 x 
10~ 16 J • m~ 3 • K -4 is the blackbody constant. 

We deduce from the constancy of sa 3 that g s ( aT ) 3 remains constant. Hence the temperature T of 
the universe behaves as <i 1 whenever g s remains constant. We conclude that as we go back in time (a 
decreases), temperature increases, as well as energy density. The early universe is hot and dense. It might 
even reach a stage where our equations no longer apply because it becomes infinitely hot and dense: this 
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is the initial singularity, sarcastically (but successfully) called big bang by Fred Hoyle, on a BBC radio 
show in 1949. 

We note that some caution has to be paid whenever some species drop out of thermal equilibrium. 
Indeed, a given species drops out of equilibrium when its interaction rate T drops below the expansion 
rate H. For example neutrinos decouple at temperatures below 1 MeV. Their temperature continues to 
decrease as a -1 and thus remains equal to T. However when kT drops below 2 m e , electrons annihilate 
against positrons with no possibility of being regenerated and the entropy of the electron-positron pair's 
is transferred to the photons. Since <y s |_ e ± = 2 + 4 • 7/8 = 11/2 and g s | = 2, the temperature of the 
photons becomes multiplied by a factor (11/4) 1 / 3 . Since the neutrinos have already decoupled, they are 
not affected by this entropy release and their temperature remains untouched. Thus we have 


T 



(51) 


We can then compute the value of g s for temperatures much smaller than m e : g s = 2 + (7/8)6(4/11) = 
3.91. We deduce that, at present time (To = 2.725 K), so/k = 2890 crrU 3 . 


As we have seen in (30), it follows from the Friedmann equation that, if the Universe is dominated 
by a component of equation of state p = wp, then the cosmic scale factor a(t) varies with time as 
t 2 /[ 3 ( 1+1 ")]. We start at time to with the energy budget: U M = 0.3, Ha = 0.7, ~ 0 (see next 

Chapters). Radiation consists of photons and relativistic neutrinos. Since generically 


Pr 


where g(T) is the effective number of degrees of freedom 

s(t ) = V * 


E 

bosons i 


I) 


+UL 


9i 


I) 


\ 4 


(52) 


(53) 


(we have taken into account the possibility that the species i may have a thermal distribution at a temper¬ 
ature Tj different from the temperature T of the photons), we have 


P R (t) = (hi*) 




(54) 


where N* el (t) is the number of relativistic neutrinos at time t. At present time to, we have H 7 = 
Pi (to)/Pc = 2.48 x 10 -5 hr 2 and the mass limits on neutrinos imply l\T el (to) < 1. In any case, 

n R « n M . 


For redshifts larger than 1, the cosmological constant becomes subdominant and the universe is 
matter-dominated (a(t) ~ t 2,/;i ). In the early phase of this matter-dominated epoch, the Universe is 
a ionized plasma with electrons and protons: it is opaque to photons. But, at a time i rec , electrons 
recombine with the protons to form atoms of hydrogen and, because hydrogen is neutral, this induces the 
decoupling of matter and photon: from then on (i rec < t < to), the universe is transparent] This is the 
important recombination stage. After decoupling the energy density p 7 ~ T 4 of the primordial photons 
is redshifted according to the law 


m 

To 


_ao_ 

a(t ) 


1 + 2 


(55) 


6 After recombination, the intergalactic medium remains neutral during a period often called the dark ages, until the first 
stars ignite and the first quasars are formed. The ultraviolet photons produced by these sources progressively then re-ionize the 
universe. This period, called the re-ionization period, may be long since only small volumes around the first galaxies start to be 
ionized until these volumes coalesce to re-ionize the full intergalactic medium. But, in any case, the universe is then sufficiently 
dilute to prevent recoupling. 
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Table 1: The different stages of the cosmological evolution in the standard scenario, given in terms of time t 
since the big bang singularity, the energy kT of the background photons and the redshift z. The double line 
following nucleosynthesis indicates the part of the evolution which has been tested through observation. The 
values (ho = 0.7, Cl M = 0.3, = 0.7) are adopted to compute explicit values. 


t 

/cT 7 (eV) 

z 


c 

O 

l 

1 — 1 

Oi 

O 

>-$ 

2.35 x 10“ 4 

0 

now 

~ Gyr 

~ 10“ 3 

~ 10 

formation of galaxies 

~ 4 x 10 5 yr 

0.26 

1100 

recombination 

t eq ~ 4 x 10 4 yr 

0.83 

3500 

matter-radiation equality 

3 min 

6 x 10 4 

2 x 10 8 

nucleosynthesis 

1 s 

10 6 

3 x 10 9 

e + e _ annihilation 

4 x 10 -6 s 

4 x 10 8 

10 12 

QCD phase transition 

< 4 x 10 -6 s 

> 10 9 


baryogenesis 




inflation 

t = 0 


oo 

big-bang 


One observes presently this cosmic microwave background (CMB) as a radiation with a black-body 
spectrum at temperature To = 2.725 K or energy kT$ = 2.35 x 10 -4 eV. 

Since the binding energy of the ground state of atomic hydrogen is Ej, = 13.6 eV, one may expect 
that the energy kT rec is of the same order. It is substantially smaller because of the smallness of the ratio 
of baryons to photons rj = rib jru, ~ 5 x 10 -10 . Indeed, according to the Saha equation, the fraction x 
of ionized atoms is given by 


n p n e _ x 2 _ 4.05c 3 / m e \ : V 2 _ Eb/kT 

nHti'y (1 — x)^ ir 2 \2irkT) 


(56) 
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Hence, because t) <C 1, the ionized fraction x becomes negligible only for energies much smaller than 
E b . A careful treatment gives kT iec ~ 0.26 eV. 

As we proceed back in time, radiation energy density increases more rapidly (as a(t) -4 ) than 
matter (a(i)~ 3 ) (as a{t) decreases). At time t eq , there is equality. This corresponds to 

1 = q(teq) = 1.68 = 4.17 X 10~ 5 

1 + Z e q &0 


where we have assumed 3 relativistic neutrinos at this time. 


As we go further back in time, we presumably reach a period where matter overcame antimatter. 
It is indeed a great puzzle of our Universe to observe so little antimatter, when our microscopic theories 
treat matter and antimatter on equal footing. More quantitatively, one has to explain the following very 
small number: 



n b - n~ b 


~ 6 x 10“ 10 , 


(58) 


n 7 re 7 

where ng = n b — n b (resp. n 7 ) is the baryon (resp. photon) number density, based on baryon b and 
antibaryon b counts. The actual number comes from the latest Planck data [ 16J. 


Sakharov [17 ] gave in 1967 the necessary ingredients to generate an asymmetry between matter 
and antimatter: 


- a process that destroys baryon number, 

- a violation of the symmetry between matter and antimatter (the so-called charge conjugation), as 
well as a violation of the time reversal symmetry, 

- an absence of thermal equilibrium. 


The Standard Model ensures the second set of conditions (CP violation which was discovered by Cronin 
and Fitch 1181 is accounted for by the phase of the CKM matrix). The expanding early universe provides 
the third condition. It remains to find a process that destroys baryon number. Different roads were 
followed: non-perturbative processes (sphaleron) at the electroweak phase transition; proton decay in the 
context of grand unified theories; decay of heavy neutrinos which leads to lepton number violation, and 
consequently to baryon number violation (leptogenesis). 


2 The days where cosmology became a quantitative science: cosmic microwave back¬ 
ground 

We have recalled briefly in Section |~i~6| the history of the Universe (see Table[l]). The very early universe 
is a ionized plasma, and thus is opaque to light. But we have seen that, soon after matter-radiation 
equality, electrons recombined with the protons to form neutral atoms of hydrogen, which induces the 
decoupling of matter and photon. From this epoch on, the universe becomes transparent to light. The 
primordial gas of photons produced at this epoch cools down as the universe expands (following [55]) and 
forms nowadays the cosmic microwave background (CMB). 

Bell Labs radio astronomers Arno Penzias and Robert Wilson were using a large horn antenna in 
1964 and 1965 to map signals from the Milky Way, when they accidently discovered the CMB [3j. The 
discovery of this radiation was a major confirmation of the hot big bang model: its homogeneity and 
isotropy was a signature of its cosmological origin. Flowever, the degree of homogeneity and isotropy of 
this radiation was difficult to reconcile with the history of the Universe as understood from the standard 
big bang theory: radiation coming from regions of the sky which were not supposed to be causally 
connected in the past had exactly the same properties. 

It was for such reasons that the scenario of inflation was proposed [4j5j: an exponential expansion 
of the Universe right after the big bang. In the solution proposed by A. Guth (5], the set up was the 
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Fig. 4: Blackbody spectrum of the CMB as observed by the instrument FIRAS onboard the COBE satellite [21] 


spontaneous breaking of the grand unified theory: the corresponding phase transition was providing the 
vacuum energy necessary to initiate such an exponential expansion. This scenario proved to be difficult 
to realize and it was followed by many variants: new inflation 1191, chaotic inflation |20|, ... 


We have said that, for t < t iec , i.e. z > 1100, the Universe is a ionized plasma, opaque to 
electromagnetic radiation. This means that, when we observe the early Universe, we hit a “wall” at the 
corresponding redshift: the earlier Universe appears to our observation as a blackbody, and should thus 
radiate according to the predictions of Planck. This is certainly the largest blackbody that one could think 
of. This blackbody spectrum of the CMB was indeed observed j2TJ by the FIRAS instrument onboard 
the COsmic Background Explorer (COBE) satellite launched in 1989 by NASA (see Fig. [4]). 


More precisely, the CMB has the nearly perfect thermal spectrum of a black body at temperature 
T 7 = 2.725 K (corresponding to a number density n 7 = 411 cm -3 ): 


dp 1 = 


2 hf 


1 


4:17 f 2 df 


e hf/kT\ _ i c 3 


(59) 


(the first factor accounts for the two polarizations) or 


= 3.8 x 10~ 15 J/m 3 f y- 

\/7 


e — 1 
e /// 7 _ 1 


(60) 


din/ 

where / 7 = kT^/h = 5.7 x 10 10 Hz. 

Another expectation for COBE was the presence of fluctuations in the CMB. Indeed, if inflation 
was to explain the puzzle of isotropy and homogeneity of the CMB through the whole sky, one expected 
that quantum fluctuations produced during the inflation phase would show up to some degree as tiny 
fluctuations of temperature in the CMB. Such fluctuations were discovered by the instrument DMR 
(Differential Microwave Radiometers) onboard COBE, at the level of one part to 10 5 [221 (see Fig. [ 5 ]). 

It is primarily homogeneous and isotropic but includes fluctuations at a level of 10 -5 , which are 
of much interest since they are . 


Since the days of COBE, there has been an extensive study of the CMB fluctuations to identify 
the imprints of the recombination and earlier epochs. This uses the results of ground, balloon or space 
missions (WMAP in the US and most recently Planck in Europe). We will review this in some details in 
the next Section. 
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Fig. 5: Temperature fluctuations in the CMB observed by COBE {22j (the galaxy is removed). The galaxy as well 
as the dipole component associated to the Doppler effect due to the motion of the Earth have been removed. Blue 
areas are colder than average, red areas are warmer but no 7° region varies from the mean by more than 200 /iK 
(AT/T ~ 8.10 -5 ). 


2.1 CMB 


Before discussing the spectrum of CMB fluctuations, we introduce the important notion of a particle 
horizon in cosmology. 


Because of the speed of light, a photon which is emitted at the big bang (t = 0) will have travelled 
a finite distance at time t. The proper distance dCT} measured at time t is simply given by the integral: 


/' ' cdf r 

dph(t ) = a(t) / —^ (61) 

f°° _ dz _ 

l + [fl M (l + ^) 3 + ^(l + z) 4 + fl fc (l + z) 2 + fl A ] 1/2 ’ 


where, in the second line, we have used ( |C.2[ ). This is the maximal distance that a photon (or any particle) 
could have travelled at time t since the big bang. In other words, it is possible to receive signals at a time 
t only from comoving particles within a sphere of radius di, (t). This distance is known as the particle 
horizon at time t. 


A quantity of relevance for our discussion of CMB fluctuations is the horizon at the time of the 
recombination i.e. z iec ~ 1100. We note that the integral on the second line of ( [61] ) is dominated by the 
lowest values of z: z ~ z iec where the universe is still matter dominated. Hence 

2 fir, 

) ~ s w^~ a3Mpc ' <62) 

M rec 

One may introduce also the Hubble radius 

R H (t) = H~\z) , (63) 
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which will play an important role in the following discussion. This scale characterizes the curvature of 
spacetime at the time t {see for example Exercise B.l.b). We note that the particle horizon is simply twice 
the Hubble radius at recombination, as can be checked from ( |43j ): 

~ ■ 

iV1 rec 

This radius is seen from an observer at present time under an angle 


c ) = 


R H(t lec ) 
dA(t iec ) 


(65) 


where the angular distance has been defined in ( C.7| ). We can compute analytically this angular distance 
under the assumption that the universe is matter dominated {see Exercise C-l). Using (C.10|), we have 


dA(t ie c ) = 


a 0 r 


2 ^h 0 

^iec 


1 + Z re 

Thus, since, in our approximation, the total energy density Q 7 , is given by Q v; , 

0 H (t T J ~ Ui /2 /( 2z x J 2 ) ~ 0.015 rad SI 1 ll ' 2 ~ 1° Q}' 2 . 


( 66 ) 


(67) 


We have written in the latter equation fl T instead of Cl M because numerical computations show that, in 
case where Ha is non-negligible, the angle depends on + Q,\ = 0 T . 

We can now discuss the evolution of photon temperature fluctuations. For simplicity, we will 
assume a flat primordial spectrum of fluctuations: this leads to predictions in good agreement with 
experiment; moreover, as we will see in the next Section, it is naturally explained in the context of 
inflation scenarios. 

Before decoupling, the photons are tightly coupled with the baryons through Thomson scattering. 
In a gravitational potential well, gravity tends to pull this baryon-photon fluid down the well whereas 
radiation pressure tends to push it out. Thus, the fluid undergoes a series of acoustic oscillations. These 
oscillations can obviously only proceed if they are compatible with causality i.e. if the corresponding 
wavelength is smaller than the horizon scale or the Hubble radius: A = 2ira{t)/k < Rnif) or 

k > 2 '/S) ~ 4-1/3 ' <68) 


Starting with a flat primordial spectrum, we see that the first oscillation peak corresponds to A ~ 
Rjj( free ), followed by other compression peaks at Ru(t THr )/n (see Fig. [6]). They correspond to an 
angular scale on the sky: 


RH(t rec ) 1 = 0H{t lec ) 
dA(t Tec ) n n 


Since photons decouple at t r 
temperature^ 


we observe the same spectrum presently (up to a redshift in the photon 


Experiments usually measure the temperature difference of photons received by two antennas 
separated by an angle 9, averaged over a large fraction of the sky. Defining the correlation function 


/ AT AT \ 

cw = W (ni) ^ (n2) ) 


(70) 


1 A more careful analysis indicates the presence of Doppler effects besides the gravitational effects that we have taken into 
account here. Such Doppler effects turn out to be non-leading for odd values of ri¬ 


ll 

















Fig. 6: Evolution of the photon temperature fluctuations before the recombination. This diagram illustrates that 
oscillations start once the corresponding Fourier mode enters the Hubble radius (these oscillations are fluctuations 
in temperature, along a vertical axis orthogonal to the two axes that are drawn on the figure). 


averaged over all ni and ri 2 satisfying the condition ni • ri 2 = cos 6, we have indeed 


T(ni) — T(n 2 )' 


= 2 (67(0)-67(0)). 


(71) 


We may decompose 67(0) over Legendre polynomials: 

OO 

C(0) = —J2(2l + l)CiPi(cos9). 


(72) 


The monopole (l = 0), related to the overall temperature To, and the dipole (l = 1), due to the Solar 
system peculiar velocity, bring no information on the primordial fluctuations. A given coefficient 67; 
characterizes the contribution of the multipole component l to the correlation function. If 9 <C 1, the 
main contribution to 67; corresponds to an angular scaled ~ ir/l ~ 200° /l. The previous discussion 
(see (67) and (69)) implies that we expect the first acoustic peak at a value l ~ 200 


The power spectrum obtained by the Planck experiment is shown in Fig. [7] One finds the first 
acoustic peak at l ~ 200, which constrains the ACDM model used to perform the fit to 0 7 , = Q M + Q \ ~ 


1. Many other constraints may be inferred from a detailed study of the power spectrum [16 231. 


2.2 Baryon acoustic oscillations 

We noted in the previous section that, before decoupling, baryons and photons were tightly coupled and 
the baryon-photon fluid underwent a series of acoustic oscillations, which have left imprints in the CMB: 
the characteristic distance scale is the sound horizon, which is the comoving distance that sound waves 
could travel from the Big Bang until recombination at z = z*: 


r s 



+ 



c s (z) 

H(z) 


dz , 


(73) 


8 The Ci are related to the coefficients a; m in the expansion of A T/T in terms of the spherical harmonics Yj m : C; = 
(|a!m| 2 )m- The relation between the value of l and the angle comes from the observation that Yi m has (l — m) zeros for 
— 1 < cos 9 < 1 and ReO)™) m zeros for 0 < (j> < 2ir. 
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Fig. 7: This figure compares the temperature spectrum for the best fit ACDM model, in red, with the temperature 
angular power spectrum observed by the Planck collaboration (in blue, averaged over bins of width A l ~ 31 with 
1 a errors). The gray dots are the unbinned data. In the lower panel, the green lines show the ±ler errors on the 
individual power spectrum estimates. HD 


where c s is the sound velocity. This distance has been recently measured with precision by the Planck 
collaboration to be r s = 144.96 ± 0.66 Mpc 1161. 

We have until now followed the fate of photons after decoupling. Similarly, once baryons decou¬ 
ple from the radiation, their oscillations freeze in, which leads to specific imprints in the galaxy power 
spectrum, such as the characteristic scale r s . Indeed, remember that, until recombination, baryons and 
photons were tightly coupled (but not dark matter). Thus a given matter density perturbation may have 
travelled a distance r s in the case of baryons under the influence of radiation pressure (to which the 
photons are sensitive), whereas it did not move in the case of dark matter. This will lead, once (dark and 
baryonic) matter has collapsed into galaxies, to a secondary peak a distance r s away in the distribution 
of separations of pairs of galaxies. In 2005, Eisenstein and collaborators 1241, using data from the Sloan 
Digital Sky Survey, have indeed identified such a baryon acoustic peak in the matter power spectrum 
(Fourier transform of the two-point correlation function) on scales of order 105 h~ l ~ 150 Mpc. Fig¬ 
ure 41 shows how the ACDM model just described fares with respect to observations by comparison 


with a model with = 0, a small value of the Hubble parameter and a small fraction of matter which 
does not cluster on small scale (relic neutrinos or quintessence). 


The acoustic peak provides a standard ruler which can be used for measuring distances: measure¬ 
ments along the line of sight depend on H(z)r s whereas measurements transverse to the line of sight 
depend on the angular diameter distance d,A{z)/r s . In fact, because the analysis rests on spherically- 
averaged two-point statistics, the distance scale determined is dy defined as: 


/ rfy(z) \ 3 = / d A (z) \ 2 cz 
\ r s ) ~ \ r s J H(z)r s 


(74) 


Measuring the acoustic scale at z = 0 provides a standard distance ruler which allows to identify /7 () , 
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Fig. 8: Correlation function in redshift space for the best-fit power-law ACDM model 
best-fit model with (>\ = 0, Hq = 46 km/s/Mpc and a relic neutrino component 
quintessence component ( .1q = 0.12 (dot-dashed line) 1251. 


(dotted line) and for the 
= 0.12 (solid line) or a 


and then Q M (from U A/ lr). Going to higher z, this allows to put constraints on the recent history of the 
Universe, and thus on the evolution of the dark matter component. 


2.3 Inflation 

The inflation scenario has been proposed to solve a certain number of problems faced by the cosmology 
of the early universe [5j. Among these one may cite: 


- the flatness problem 

If the total energy density pr of the universe is presently close to the critical density, it should have 
been even more so in the primordial universe. Indeed, we can write (25 1 of Section[2]as 


Pr(t) _ i = A 
p c (t) a 2 


(75) 


where p c (t) = 3H 2 (t)/(8irG N ) and the total energy density pr includes the vacuum energy. If 
we take for example the radiation-dominated era where a(t) ~ t 1//2 , then (75 I can be written as 

(a ~ t -1 / 2 ~ a -1 ) 


PT(t) 1 _ Prjtjj) 1 

Pcit) [ p c (t v ) 



' ptK) 1 fVT\ 2 

. Pc(t v ) J V kT J 


(76) 


where we have used the fact that T{t) oc a(i ) -1 and we have taken as a reference point the 
epoch t v of the grand unification phase transition. This means that, if the total energy density 
is close to the critical density at matter-radiation equality (as can be inferred from the present 
value), it must be even more so at the time of the grand unification phase transition: by a factor 
(leV/10 16 GeV)" ~ 10~ 50 ! Obviously, the choice k = 0 in the spatial metric ensures p T = p c but 
the previous estimate shows that this corresponds to intial conditions which are highly fine tuned. 
- the horizon problem 
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We have stressed in the previous Sections the isotropy and homogeneity of the cosmic microwave 
background and identified its primordial origin. It remains that the horizon at recombination is 
seen on the present sky under an angle of 2°. This means that two points opposite on the sky were 
separated by about 100 horizons at the time of recombination, and thus not causally connected. It is 
then extremely difficult to understand why the cosmic microwave background should be isotropic 
and homogeneous over the whole sky. 

- the monopole problem 

Monopoles occur whenever a simple gauge group is broken to a group with a 17(1) factor. This is 
precisely what happens in grand unified theories. In this case their mass is of order M v / g 2 where 
g is the value of the coupling at grand unification. Because we are dealing with stable particles 
with a superheavy mass, there is a danger to overclose the universe, i.e. to have an energy density 
much larger than the critical density.. We then need some mechanism to dilute the relic density of 
monopoles. 


Inflation provides a remarkably simple solution to these problems: it consists in a period of the 
evolution of the universe where the expansion is exponential. Indeed, if the energy density of the universe 
is dominated by the vacuum energy p vac (or by some constant form of energy), then the Friedmann 
equation reads 



PvSLC 

3 m 2 p 


(77) 


where m p = (8 ttG n ) 1//2 is the reduced Planck mass. If p vac > 0,this is readily solved as 


a(t) = 77, 


-i „Hv, 
vac 0 


with H v 



(78) 


Such a behaviour is in fact observed whenever the magnitude of the Hubble parameter changes slowly 
with time i.e. is such that H | <C H 2 . 

As we have seen in (jTOj), such a space was first proposed by de Sitter fl0||TT[ with very different 
motivations and is thus called de Sitter space. 

Obviously a period of inflation will ease the horizon problem. Indeed, the particle horizon size 
during inflation reads, following ( [61] ) 

, ... ,,, f l cdt! 

dph{t) Ide Sitter = a W 


ti a O') 


TTv 


e H vac (t u) for H ^ t _ > L 


(79) 


It follows that a period of inflation extending from ti to tf = U + At contributes to the particle horizon 
size a value ce HvacAt /H vac , which can be very largcj^] 

We note that de Sitter space also has a finite event horizon. This is the maximal distance that 
comoving particles can travel between the time t where they are produced and t = oo (compare with 

m 

cdt' 

d eh {t) = a{t) / — . (80) 

Jt a \t') 

In the case of de Sitter space, this is simply 


c 

dehit) | de Sitter = 77 = , (81) 

vac 

i.e. it corresponds to the Hubble radius (constant for de Sitter spacetime). This allows to make an 
analogy between de Sitter spacetime and a black hole: we will see in Section [3T2] that a Schwarzschild 

9 We also note that, in a pure de Sitter space, the particle horizon diverges as we take ti —> —oo. This reflects the fact that, 
in a de Sitter space, all points were in causal contact. 
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black hole of mass M has an event horizon at the Schwarzschild radius Rs = 2G N M (see also Exercise 
1-3 of Section[l]for a comparison between the Schwarzschild and the de Sitter metric in its static form). 
Thus, just as black holes evaporate by emitting radiation at Hawking temperature Th = l/(4irRs) (see 
Eq. ( [ 1 36| >), an observer in de Sitter spacetime feels a thermal bath at temperature Th = H/(2tt). 

We see that it is the event horizon that fixes here the cut-off scale of microphysics. Since it is equal 
here to Hubble radius, and since the Hubble radius is of the order of the particle horizon for matter or 
radiation-dominated univers^J it has become customary to compare the comoving scale associated to 
physical processes with the Hubble radius (we already did so in our discussion of acoustic peaks in CMB 
spectrum; see Fig. [6], and Fig.[9]below). 

A period of exponential expansion of the universe may also solve the monopole problem by dilut¬ 
ing the concentration of monopoles by a very large factor. It also dilutes any kind of matter. Indeed, a 
sufficiently long period of inflation “empties” the universe. However matter and radiation may be pro¬ 
duced at the end of inflation by converting the energy stored in the vacuum. This conversion is known 
as reheating (because the temperature of the matter present in the initial stage of inflation behaves as 
a ~ 1 (t) oc e~ Hvact , it is very cold at the end of inflation; the new matter produced is hotter). If the reheat¬ 
ing temperature is lower than the scale of grand unification, monopoles arc not thermally produced and 
remain very scarce. 

Finally, it is not surprising that the universe comes out very flat after a period of exponential 
inflation. Indeed, the spatial curvature term in the Friedmann equation is then damped by a factor a 2 oc 
e - 2 H va , c At p or cxani p| e ^ a value H vac At ~ 60 (one refers to it as 60 e-foldings) would easily account 
for the huge factor 10 50 of adjustment that we found earlier. 


Most inflation models rely on the dynamics of a scalar field in its potential. Inflation occurs 
whenever the scalar field evolves slowly enough in a region where the potential energy is large. The set 
up necessary to realize this situation has evolved with time: from the initial proposition of Guth |5J where 
the field was trapped in a local minimum to “new inflation” with a plateau in the scalar potential |19 261, 
chaotic inflation |20| where the field is trapped at values much larger than the Planck scale and more 
recently hybrid inflation [27] with at least two scalar fields, one allowing an easy exit from the inflation 
period. 

The equation of motion of a homogeneous scalar field <p(t) with potential V{<j>) evolving in a 
Friedmann-Robertson-Walker universe is: 


4> + 3H<j) = (82) 

where V'{(j>) = dV/dcj). The term 3 Hcf) is a friction term due to the expansion. The corresponding energy 
density and pressure are: 


P = + 

v = ^ 2 - h (</>). 


(83) 

(84) 


We may note that the equation of conservation of energy p = —3 H(p + p) takes here simply the form 
of the equation of motion ([82]). These equations should be complemented with the Friedmann equation 

m- 

When the field is slowly moving in its potential, the friction term dominates over the acceleration 
term in the equation of motion (82) which reads: 

3 Hj>~-V'(<l>) . (85) 

10 In an open or flat universe, the event horizon 80 is infinite. 
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(86) 


The curvature term may then be neglected in the Friedmann equation (|77]) which gives 


H 2 ~ 


P 


V 


3m? 3m? 


Then the equation of conservation p = —3H(p + p) = —3 Hcj? simply gives 

rl? 

H ~ - 


2 mi 


(87) 


It is easy to see that the condition |iT| <C H 2 amounts to 4?/2 <C p/3 ~ V{4>)/3, i.e. akinetic energy 
for the scalar field much smaller than its potential energy. Using ( [85] ) and ( |86| ), the latter condition then 
reads 


e 


I (m p V' 
V 


f\ 2 


< 1 


( 88 ) 


The so-called slow roll regime is characterized by the two equations (85 l and (86), as well as the condition 
([88]). It is customary to introduce another small parameter: 


m 2 V" 

^ ^ v « 1 


(89) 


which is easily seen to be a consequence of the previous equation ^ 1 2 [ 

An important quantity to be determined is the number of Flubble times elapsed during inflation. 
From some arbitrary time t to the time t e marking the end of inflation (i.e. of the slow roll regime), this 
number is given by 


N(t) = / H{t)dt . 


(92) 


It gives the number of e-foldings undergone by the scale factor a(t) during this period (see (78). Since 
dN = —Hdt = — Hd(j)/4’ one obtains from (85 l and (86) 


ft IV 

N (4>) = / —TF, d $ ■ 
J6 r m? V' 


(93) 


During the inflationary phase, the scalar fluctuations of the metric may be written in a conformal 
Newtonian coordinate system as: 


ds 2 = a 2 [(1 + 2 <&)dr] 2 - (1 - 2<S>)6 ij dx i dx j ] 


(94) 


where rj is conformal time ( adr/ = dt = da/a). We may write the correlation function in Fourier space 
V s {k )by 

= 2t T 2 k- 3 V s (k)5 3 (k - k 7 ) . (95) 

The origin of fluctuations is found in the quantum fluctuations of the scalar field during the de Sitter 
phase. Indeed, if we follow a given comoving scale a(t)/k with time (see Fig. [9]), we have seen in 
Section 2.1 that, some time during the matter-dominated phase, it enters the Hubble radius. Since a(t) 


‘Differentiating 1 85 one obtains 


V = ■ 


(90) 


12 


Note that one finds also in the literature the slow roll coefficients defined from the Hubble parameter [281 

2 


e H = 2m, 


( H'W 


_ 0 2 H'\4>) 

= e , Vh = 2m p -=r 1 -e. 


(91) 
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Fig. 9: Evolution of a physical comoving fluctuation scale with respect to the Hubble radius during the inflation 
phase {Ru{t) = H~. a J.), the radiation dominated phase (Rn(t) = 2 1) and matter dominated phases (Rn(t) = 
3t/2). 


is growing (even exponentially during inflation) whereas the Hubble radius is constant during inflation, 
this means that at a much earlier time, it has emerged from the Hubble radius of the de Sitter phase. In 
this scenario, the origin of the fluctuations is thus found in the heart of the de Sitter event horizon: using 
quantum field theory in curved space, one may compute the amplitude of the quantum fluctuations of 
the scalar field; their wavelengths evolve as a(t)/k until they outgrow the event horizon i.e. the Hubble 
radius; they freeze out and continue to evolve classically. The fluctuation spectrum produced is given by 


V s (k) = 


H 2 

A& 


H 

2tt 


127r 2 m® 

k=aH p 


V :i 

yj2 


(96) 


k=aH 


where the subscript k = aH means that the quantities are evaluated at Hubble radius crossing, as ex¬ 
pected. We also note that H i.e. Rh sets the scale of quantum fluctuations in the de Sitter phase (see 
Appendix\D[for details, in particular Eq. (DAO )): Rh is indeed the dynamical scale associated with the 
physics of fluctuations (which happens to coincide with either of the kinematical scales which are the 
particle horizon for the matter dominated phase, or the event horizon in the inflation phase). 

The scalar spectral index ns(k) is computed to be (see Exercise D-l in Appendix |d|): 


ns(k) - 1 = 


dlnVs(k) 

dink 


= — 6 e + 2?7 . 


(97) 


Thus, because of the slow roll, the fluctuation spectrum is almost scale invariant, a result that we have 
alluded to when we discussed the origin of CMB fluctuations. One of the highlights of the Planck 
cosmology results |16]29j is the confirmation that the spectrum is not scale invariant i.e. ns is different 
from 1: 

n s = 0.9603 ± 0.0073 . (98) 

In other words, we are really in a slow roll phase, i.e. an unstable phase which is crucial since eventually 
one has to get out of inflation and reheat. 

The observation by the COBE satellite of the largest scales has set an important constraint on 
inflationary models by putting an important constraint on the size of fluctuations (see the caption of 
Figure [5]). Specifically, in terms of the value of the scalar potential at horizon crossing, this constraint 
known as COBE normalization, reads (see (|96|): 


1 ys/2 


777,3 yl 


= 5.3 x 1(T 4 . 


(99) 
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Using the slow roll parameter introduced above in ( 881 , the COBE normalization condition can be written 


as 


U 1/4 ~ 0.03 e 1/4 m p . 


( 100 ) 


Besides scalar fluctuations, inflation produces fluctuations which have a tensor structure, i.e. pri¬ 
mordial gravitational waves. They can be written as perturbations of the metric of the form 


ds 2 = a 2 [q lw + hjj] , (101) 

where is a traceless transverse tensor (which has two physical degrees of freedom i.e. two polariza¬ 
tions). The corresponding tensor spectrum is given by 


V T {k) 



with a corresponding spectral index 


n T (k) = 


dluVrik) 
d In k 


= —2e . 


( 102 ) 


(103) 


We note that the ratio Vt/Vs depends only on oi 2 /// 2 and thus on e, which yields the consistency 
condition: 



8<j) 2 

m 2 p H 2 


= 16e = —8 riT ■ 


(104) 


2.4 Inflation scenarios 

We conclude this discussion by reviewing briefly the main classes of inflation models. Let us note 
that, for an inflationary model, the whole observable universe should be within the Hubble radius at the 
beginning of inflation. This corresponds to a scale 


k = a 0 H 0 = aH\ hc 


(105) 


where h.c. stands for “horizon crossing”. This puts a constraint on the number of e-foldings ( [93] ) between 
horizon crossing and the end of inflation (i.e. end of the slow roll regime) necessary for the inflation to 
be efficient. More generally, one defines [30] N(k) as the number of e-foldings between the time of 
horizon crossing of the scale k ( tk,h.c . at which k = aH(tk,h.c.)) an< i the end of inflation (f e ): 


N(k) = In (a(t e )/a(tk,h.c.)) 


(106) 


Distinguishing the time when the universe reheats (t r h) and the time of matter-radiation equality ( t eq ), 
we have 


®(^fc,ti.c.) _ g— N(k) tt(^e) 0>(trh) ®(^eq) 

(30 O’itrh) O’iteq) 

_ e ~N(k) f P(krhY\ ^ f P(teq) \ ^ a (t e g) 

V P(te) ) \P(trh)J “0 

C -N(k) f P(Uh) \ 1/3 f P0 \ V4 ( a(teq ) \ V4 
\p(te)J \p(trh)J \ (t0 / 


(107) 


where we used p(t eq ) = po(ao/a(t eq )) 3 and we have assumed matter domination between the end of 
inflation and reheating. Using k = aH(tk,h.c.) an d H{tk,h.c.)/Ho = (pk,h.c./ Po) 1 ^ 2 - one obtains from 
(|106|), with transparent notations: 


N(k) = 62 - In 


a^Ho 


, 10 16 GeV , 

— In- — -h In 


V, 


1/4 


Kh.c. 1 Ve /A 

-In 


k,h.c. 


V e 


1/4 


1/4 

Prh 


(108) 
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Fig. 10: Regions corresponding to the different inflation models in the plot r vs n s HD 


In the rather standard case where Vkji.c. ~ V e ~ p r h ~ (lO 16 GeV) 4 , this requires 60 e-folding for 
inflation to be efficient at the scale of the observable Universe. A length scale corresponding to 200 Mpc 
(hence k = 27t/200 Mpc) corresponds to N(k) ~ 50. 

The three main classes of inflation models (see Figs. [T0| and [TT]) are: 

- convex potentials or large field models (0 < rj < 2e) 

The potential is typically a single monomial potential: 

= , (109) 

with n > 1. Since the slow roll parameters are e = n 2 (m p /(j)) 2 /2 ancj^j// = n(n — l){m p / cj)) 2 , 
the slow roll regime corresponds to cj) 3> m p \Jn(n - 1) (for n > 2). Because the field has 
a value larger than the Planck scale (hence the name “large field model”), this might seem out 
of the reach of the effective low energy gravitation theory. But A. Linde |20| argued that the 
criterion is rather V < m p ; in fact, he suggested that the scalar field emerges from the Planck 
era with a value fa such that V(fa) ~ m A p i.e. fa ~ m p / A 1 /”. This corresponds to the chaotic 
inflation scenario, the simplest example of which being a quadratic potential J20| . A difficulty 
is that the COBE normalisation imposes an unnaturally small value for the A coupling: A ~ 
(5.3 x 10 _4 n) 2n A ri_2 ^. Another drawback is the large value of the field which makes it necessary 
to include all non renormalisable collections of order (0/M p ) n+p , unless they are forbidden by 
some symmetry. 

The limit case in this class is the exponential potential 

U(</>) = Vbexp(—A 4>/m p ) (110) 

which leads to a power law inflation [32]: a(t ) oc t 2 ^ 2 . This model yields rj = 2e = A 2 , hence 
r = —8 (ns — 1). It is incomplete since inflation does not end. 

- concave potential or small field models (r\ < 0) 

In this class, illustrated first by the new inflation scenario ]19[|26| , the field 4> starts at a small value 
and rolls along an almost flat plateau (where V"(cj)) < 0) before falling to its ground state. This 
type of potential, often encountered in symmetry breaking transitions may be parametrized, during 

13 Note that e < r/ < 2e. 
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Fig. 11: Constraints set by Planck data on various inflation models in the plot r (evaluated at the pivot scale 
fc* = 0.002 Mpc -1 ) vs n s |29). Small dots correspond to models with 50 e-foldings, large dots to 60 e-foldings. 


the phase transition by: 

v{4>) = v 0 

where the dots indicate higher order terms not relevant for inflation. In the same class appears the 
so-called “natural inflation” potential [ |33[|34J 

= Vb [1 + cos(0//)] . (112) 



A difficulty shared by the class of small field models is the unnaturalness of the initial conditions: 
why start at the height of the potential, in a plateau region or close to an unstable extremum? In 
the case of a symmetry breaking potential, the rationale could be thermal: the restoration of the 
symmetry at high temperature naturally leads to start at the unstable “false vacuum”. 

The first model proposed for inflation [4] was based on a modification of gravity described by the 
following action: 

f / \ 

S = ^ l d A x^j[R + aR 2 J, (113) 


where R is the Ricci scalar associated with the metric g inj . It may be proved that this is equivalent 
to standard Einstein gravity plus a scalar field with a potential that falls in the same class as we 
just discussed (see Exercise 2-1). 


- hybrid models (0 < 2e < rj) 

The field rolls down to a minimum of large vacuum energy (where V"((f>) < 0) from a small initial 
value. Inflation ends because, close to this minimum, another direction in field space takes over 
and brings the system to a minimum of vanishing energy |27[. Such models, which thus require 
several fields, were constructed in order to allow inflation at scales much smaller than the Planck 
scale. In this case (see (100 1 , e can be very small and thus r as well. This class of model is the 
only one that can accomodate values of ns larger than 1 (i.e. a blue spectrum). 

We note that, for most models, the inflation scale is much larger than the TeV scale. This should 
in principle lead us to consider inflation models in the context of supersymmetry, in order to avoid 
an undesirable fine tuning of parameters. In fact, supersymmetric (and superstring) theories are 
plagued with the presence of numerous flat directions: this might be a blessing for the search of 
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inflation potential] Two such types of potentials rely on the properties of basic supersymmetry 
multiplets: they are called F-tcrm |35|[36| and /9-term |37 381 inflation and fall in the category of 
hybrid inflation. They are not presently favoured by Planck data because they give too large values 
of ns (of the order of 0.98). 


Exercise 2-1 : We show that a certain class of 4-dimensional models which extend Einstein theory 
are equivalent to Einstein gravity with a scalar field coupled to the metric (i.e. a scalar-tensor theory). 
This class of models is described by the following action: 

yyjf r _ _ 

s = ~Y J d A xy/^f(R) + Sm(f>,g^) , ( 115 ) 

where R is the Ricci scalar associated with the m etric (j iw , and f(R) is a general function of this Ricci 
scalar. We note that the Starobinsky action (113) corresponds to f(R) = R + aR 2 . 


Let us consider the more general action: 


S = ^ d A x\/-c 


df 


/(x) + ^-x) 


T S m (ppi 9pv 


(116) 


a) Show that the variation o/( |77(5| ) with respect to x leads to (115). 

b) Redefining the metric and the scalar field through 


df 


9pu — , Su v if — \ rn p log 


d\ 


df 


dx 


show that one recovers the familiar form of the scalar-tensor gravity: 

,2 


S = J d A x\pRj 


m " i 

-fR--dyd^-V{f) 




(117) 


(118) 


where one will express the potential V(6) in terms ofx, f(x ar, d df /dx, and give the explicit form of 

c) Identify the potential V (f) in the case of the Starobinsky function f(x) = X + a X 2 - 
Hints: a) x = R> under the condition that f" 7 / 0 . 

_ ^,2 Xdf/dx-f 


b)V = m 2 p and A(f) = exp [f/(m p V 6 )]. 


3 Light does not say it all (1): the violent Universe 

The Universe is the siege of many violent phenomena; one may cite explosions like supernovae, gamma 
ray bursts (GRB) or the emission of energetic particles by active galaxy nuclei (AGN), quasars, blazars... 
The time constants t associated with the phenomena are very short on the scale of the Universe. For 
example, a GRB may be visible on the sky only for a few seconds. This means that the distance scales cr 
involved are very small: the distance that light travels in 10 seconds is only 3 million km, that is 0.002 
astronomical unit (1 a.u. is the Sun-Earth distance). Indeed, very compact objects, such as neutron stars 
or black holes, are at the heart of such violent phenomena. We will start by reviewing the origin of such 
compact astrophysical objects, which appear at the end of the life of a star. 

14 One possible difficulty arises from the condition ( |<S9| ) which may be written as a condition on the mass of the inflaton field 

m 2 <g H 2 . (114) 

Any fundamental theory with a single dimensionful scale (such as string theory) runs into the danger of having to fine tune 
parameters in order to satisfy this constraint. This is known as the p problem. 
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3.1 The end of the life of a star: from white dwarfs to neutron stars and black holes 

The evolution of a generic gravitational system such as a star is governed by two competing processes: 
gravitational forces which tend to contract the system and thermal pressure which is due to the ther¬ 
monuclear reactions within, which tendto expand the sytem. In a stable star like our Sun at present, the 
two processes balance each other. But when the nuclear fuel is exhausted, the (core of the) star starts to 
collapse under the effect of gravity; the gravitational energy thus released heats up the outer layers of the 
star, which produces the explosive phenomena that we observe. 

But what is the fate of the collapsing core? Gravitational pressure is eventually counterbalanced 
by quantum degeneracy pressure. Let us explain the nature of this pressure. Since matter is made of 
fermions of spin 1/2, Pauli principle applies: two fermions cannot be in the same state. Fermionic 
matter will thus resist at some point to excessive pressure. 

Let us be more quantitative. Since there are 4irp 2 dp/ ( 2 tt/j ) :! levels per unit volume with momentum 
between p and p + dp and two spin states per level, the number of fermions per unit volume is given in 
terms of the maximal momentum by 


n = 


2 

(27 xh) 3 


fp F 


*iirk 2 dk = 


Pf 

3ir 2 h 3 


(119) 


The energy of the highest level, or Fermi energy f jr, is therefore given in terms of the number density n. 
If the particles are non-relativistic, then 


e F = — = - (Svr 2 ) 
2m 2 [ J 




n 2/S 


m 


( 120 ) 


On the other hand, the gravitational energy per nucleon of a system of size R and mass M (with 
N = 4irR 3 riN /3 = M/m^r nucleons) is 


G N Mm N 2 N 

= G N m N^ = Hr 


(y) G N m%N 2 / 3 nft 3 . 


y R " JV R 

The Fermi energy starts to dominate over the gravitational energy for 


H 3 >___ 

N (3tT 2 ) 2 /3 


G N m 2 N m \ N 2/3 2/3 
h 2 1 ’ 


( 121 ) 


( 122 ) 


(123) 


where v = njy/n (u depends on the species of the fermions that are degenerate; see below), or 

RM l ' 3 < — — rHV 2 / 3 , 
ac me 

where, as above, a G = ( G N m 2 N /hc ) ~6 x 10 -39 . 

We see from (122) that gravitational collapse is first stopped by the quantum degeneracy of elec¬ 
trons: the corresponding astrophysical objects are known as white dwarfs. Writing thus m = m e and 
v = 2 (two nucleons per electron), we find that RM 1 / 3 10 2 RqMq 3 . A white dwarf with M = M & 

has radius R ~ 10 2 R. . and density p ~ 10 6 /?©. It is more compact than a star. 

If density contin ues to increase, the value of the Fermi energy is such that the fermions are rela¬ 
tivistic: it follows from (119) that pp > me reads (37T 2 ) 1 / 3 hn 1 ^ 3 > me or, using n = 3 N/ (47 ti 'R 3 ), 


-Y /3 


4 J 


me* 


(124) 
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But, since e F ~ p F c = (3tt“) hen 1 both e F and e g scale like n 1 /' 3 . Quantum degeneracy pressure 

can overcome gravitational collapse only for N < 3 v /^oy/ !/ ' 2 /(2z^ 2 ), or 

M < 3\frraQ 3 ^ 2 mN/(2n 2 ) ~ 1 Mq/u 2 . (125) 


This bound is the well-known Chandrasekhar limit for white dwarf masses (a more careful computation 


gives a numerical factor of 5.87 [39]). The radius of the object then satisfies (see ( 124 1) 


30r —1/2 he 1 
R < —-— a, - 


k g 


me 2 v 


(126) 


Setting m = m e gives a limit value of some 10 4 km. 

For even higher densities, most electrons and protons are converted into neutrons through inverse 
beta decay (p + e~ —yn + v). A new object called neutron star forms when the neutron Fermi energy 
balances the gravitational energy. Writing m = m n instead of m e in (123 1 , we now have RM 1 ,/;i ~ 


10 5 RqMq , c> : a neutron star with M = M 0 has radius R ~ 10 5 i? 0 and density p ~ 10 1 


r l/3. 


-5, 


ll5, 


The bound d 125 l obtained above in the case of relativistic fermions (neutrons in this case) is called 
the Oppenheimer-Volkoff bound: more precisely, the maximal mass of a neutron star is M = 0.7 M 0 , 


with a corresponding radius R = 9.6 km (cf. (1261 with rn = m n ). If the mass is larger, the star 
undergoes gravitational collapse and forms a black hole. 


3.2 Gravitational collapse: black holes 

Let us first backtrack a little and return to Einstein’s equations Q. Because they are non-linear, there 
are few solutions known. The first exact non-trivial solution was found in late 1915 by Schwarzschild, 
who was then fighting in the German army, within a month of the publication of Einstein’s theory and 
presented on his behalf by Einstein at the Prussian Academy in the first days of 1916 [9|, just before 
Schwarzschild death from a illness contracted at the front. It describes static isotropic regions of empty 
spacetime, such as the ones encountered in the exterior of a static star of mass M and radius R. 

The Schwarzschild solution reads, for r > R {see Exercise 3-1), 

l - dt 2 -(l~ ^ dr 2 _ r 2 dd 2 _ r 2 sin 2 M(j) 2 (127) 

r ) V r ) 

The Schwarzschild solution is singular at r = Rs = 2G N M, a distance known as the Schwarzschild 
radius. This is not a problem as long as Rs < R since this solution describes the exterior region of the 
star. A different metric describes the interior. On the other hand, we will see in Section [T2| that, in the 
case where R < Rs, i.e. 2 G N M/R > 1, the system undergoes gravitational collapse and turns into a 
black hole. 



Exercise 3-1 : In this exercise, we derive the Schwarzschild solution (127 1. Because we look for 
static isotropic solutions, we may always write the spacetime metric 


ds 2 = e 2u{r) dt 2 - e 2X{r) dr 2 - r 2 ( dO 2 + sin 2 Odf 2 ) . 

In other words, the only non-vanishing elements of the metric are: 

9tt = e 2u(r \ g rr = -e 2A(r) , g ee = - r 2 , g H = -r 2 sin 2 6 . 
a) Work out the Christojfel symbols from \B. /1 and the Ricci tensor components from (B.3 1. 


(128) 


( 129 ) 


5 We have absorbed a general function e 2,i T) i n front of the last term by redefining the variable r. 
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b) Show that the Einstein’s equations in the vacuum simply amount to a condition of vanishing 
Ricci tensor: 

Rfiu = 0 • (130) 

c) From the vanishing of Rtt and R rr and the fact that, at large distance from the star, space 
should be flat, both X and u should vanish at spatial infinity. Hencedecluce that 

X = —v . (131) 


d) From the vanishing of Rqq, deduce that 


2v i 

9tt = e = 1 - 


2G„M 


(132) 


Hints: a) 

Rtt = ^' + ^-AV + ^je 2 M), 

9 \' 

r> // /2 i \ / / i 

Rrr — —V — V + X V + - , 

r 

Rqq = 1 - (l + rtJ - rX 1 ) e~ 2X , 

R .00 = Ree sin 2 0 . (133) 

b) reads R IW — ^g^R = 0. Contracting with g ,iv yields R = 0. 

d) The constant of integration is identified with the mass M because, in the Newtonian limit, 
gtt = 1 + 2<I> where $ is the Newtonian potential. 


Exercise 3-2 : What is the Schwarzschild radius of the sun? of an astrophysical object of mass 
3 X 10 6 M 0 ? 


Hints: Do not forget that we have set c = 1. Otherwise, Rs = 2G N M/c 2 , that is 2.95 km for the 
sun, 8.85 x 10 9 m = 0.06 au for an object of mass 3 x 10 6 Mg, 


We now understand that, when (the core of) a star of mass M in gravitational collapse overcomes 
the degeneracy pressure of neutrons to reach a size R < Rg = 2 G N M, nothing seems to drastically 
change for observers located at distances r > Rg. However, the behaviour of the Schwarzschild met¬ 
ric [127] appears to be singular: g a vanishes and g rr diverges. It took some time (Lemaitre again!) to 


realize that this was not the sign of a real singualrity but was just an artifact of the choice of coordinates: 
other choices lead to a regular behaviour (see Exercise 3-3). The true singularity lies at r = 0 where the 
collapsing matter ends up. 


In order to understand the nature of the surface at r = Rg, let us keep for a moment longer the 
Schwarzschild coordinates and consider sending a light signal radially from some point n to > r\ 
where it is received a time At later. Since ds 2 = 0 (as well as dO = df = 0), we have simply 


h 1 ” 2 dr 

1 “ L (1 - Rs/r) ■ 


(134) 


If n < Rg, this is finite only for r 2 < Rg, in which case it is simply r 2 ~r\+Rg In [(Rg — rf)/(Rs ~ d)]- 
In other words, signals emitted from within the Schwarschild radius never reach the outside. There is 
really a breach of communication. Indeed, the surface r = Rg is an event horizon (see Section |27T|). 
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Let us take this opportunity to present a classical interpretation of the Schwarzschild radius. Re¬ 
member that the existence of black holes was conceived by Michell [ |40| and Laplace |4T| centuries 
earlier than general relativity. Indeed, the classical condition for escape a body of mass m and velocity 
v from a spherical star of mass M and radius R is 


1 2 
-mv 

2 


> 


G N mM 

R 


(135) 


Thus, not even light (v = c) can escape the attraction of the star if R < 2 G N M/c 2 , the Schwarzschild 
radius. 


We note that the Schwarzschild horizon is a fictitious surface, in the sense that an observer crossing 
this surface would not experience anything particular (we said that there exist coordinates where the 
behaviour at R$ is regular), except deformations due to tidal forces because it comes closer to a very 
massive object. But once it has crossed this fictitious surface, there is no way to backtrack: the further 
information that might be gained is lost for ever to the outside world. A useful picture is the one of a 
person swimming in a river with a waterfall downstream: swimming in the river involves no danger as 
long as one is safely far from the waterfall, but, at some point the swimmer crosses a fictitious line (the 
“horizon” of the waterfall) which is the point of no return: even the best swimmer is attracted towards 
the “singularity” of the waterfall. 


So far, our description has been purely classical. Quantum mechanical processes change this 
picture. Indeed, S. Hawking |42j pointed out that black holes emit radiation through what is known as 
the process of evaporation. Indeed, it can be shown that an accelerated observer sees a thermal bath of 
particles at a temperature T = ha/ (27r) (a being the acceleration): this is the so-called Unruh effect [431. 
Now, an observer who is at a fixed distance r > Rs from the horizon of a black hole in Schwarzschild 
coordinates has an acceleration a = Rs( 1 — Rs/R) l,/2 /(2r 2 ) (see Exercise 3-4). For r ~ Rs, it thus 
observes a thermal bath of particles at a temperature measured by an observer at infinity to be 


T h 


h 

47 tRs 


(136) 


The Hawking evaporation process is important to understand the non-observation of primoridal black 
holes, which would be due to fluctuations of density during the Planck era: such primordial black holes 
have evaporated. 


A final comment using the Schwarzschild coordinates (127): we see that, when r crosses Rs, the 
respective signs of gtt and g rr changes. In other words, t becomes a spatial coordinates and r becomes 
time: the movement towards the central singularity is the clock that ticks. 


Exercise 3-3 : Define the Kruskal coordinates (v, a. 6 , (t>) related to the Schwarzschild coordinates 
(t,r,6,(t>) through [44]: 


for r > Rs , 

u = (r/R s - 1) 1/2 e r/2Rs cosh(f/2i?s) , 
v = ( r/R s - 1) 1/2 e r/2Rs smh(t/2R s ) , 


for r < Rs , 

u = (1 — r/R s ) l/2 e r / 2Rs sinh(f/2i?s) , 



v = (1 - r/R s ) 1/2 e r/2Rs cosh(t/2R s ) . 

(137) 

Deduce from (127) the form of the metric in Kruskal coordinates: 


ds 2 = 4i? V 
r 

- r / R S ( dv 2 _ du 2 ) _ r 2 ( d0 2 + sin 2 0d ^ 

(138) 


where r is given as an implicit function of u and v: 


- 1 ) e r/Rs = u 2 - v 2 . (139) 

Rs J 
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In order to be more quantitative, let us follow the analysis of Oppenheimer and Snyder fi45\] who 
were the first to discuss the collapse into a black hole. We consider a fluid of negligible pressure, thus 
described by the energy-momentum tensor (see ( |27| )) T )W = plffifi, and study its spherically symmetric 
collapse. 


It turns out that we have already studied this system when we discussed the evolution of a homo¬ 
geneous and isotropic universe in Section 1.3 The metric is given by 


(140) 


as in ( |79|P’1 fl , /2<7 the Einstein tensor components are the same as in 
r so that a(0) = 1. Thus 

p(t) = p(0)/a 3 (t) (141) 


(B.5\B.6). We normalize the coordinate 


ds z = dt z — a z (t) 


dr 2 


1 — kr 2 


+ r 2 d& 2 + r 2 sin 2 ddf 2 


and Einstein’s equations simply read: 


a 2 + k 
cl 2 -{- 2 ci'd ~t~ k 


87 tGn p( 0) 
3 a 

0 . 


(142) 

(143) 


Assuming that the fluid is initially at rest (a = 0), we obtain from (142 1 

8 itGn 


k = 


-p{ 0)• 


Thus, (142 1 simply reads 


a 2 ' 


! (i) = k [a \t) - 1] . 
The solution is given by the parametric equation of a cycloid: 

t/j + sin fl> 


t = 


a = 


2 y/k ’ 
1 + cos 


We see that a vanishes for ip = tt, that is after a time 


it 7r 


T = 


2 Vk 2 VSttG^pCO) 


1/2 


(144) 

(145) 


(146) 


(147) 


Thus a sphere initially at rest with energy density p(0) and negligible pressure collapses to a state of 
infinite energy density in a finite time t. 


In the case of a star of radius R and mass M, this solution for the interior of the star should be 
matched with the Schwarzschild solution (127) describing the exterior. The correspondence between the 
interior and exterior coordinates is simply r = Ra(t), 6 = 6 and o = (j), with a more complicate relation 
between t and t (see Ref. /39 / section 11.9). The first relation ensures that 


= 2MG n = R s 
R 3 R 3 


(148) 


in agreement with (144) and M = (4n/3)p(0)R 3 . 


16 except that we do not normalize k to ±1 or 0 because we are looking at a different system. We will see just below that it is 
fixed by initial conditions. We add a hat to this system of coordinates to distinguish it from the Robertson-Walker coordinates, 
as well as from the Schwarzschild coordinates that we will use later. 


33 


















Exercise 3-4 : We consider an observer at rest outside the horizon of a black hole described by the 
Schwarzschild metric (127) which we write (see Exercise 3-1): 


ds 2 = e 2v( - r) dt 2 - e-^dr 2 - r 2 ( d6 2 + sin 2 ddf 2 ) . (149) 


The observer velocity is U 11 = ^ wzY/z = Sff. 

a) Show that = ^w u — wwhere w^ (the covariant derivative is defined in 

\B.4\ of Appendix^. 

b) Deduce that the acceleration A 1 ' = U f W f ,U 11 of the observer is simply A 11 = —w ,J . 

c) Show that the acceleration a 2 = —uAuif, is given for the observer at fixed r, 0 and (b by 


a = 


G n M 

r 2 (l — 2G JV M/r) 1 /2 ' 


(150) 


We finally note that, at large distance, the black hole is only caracterized by its mass M. Black 
holes arc indeed very simple objects, somewhat similar to particles: Schwarzschild black holes arc only 
characterized by their mass. Other more complex solutions were found later but it was realized that 
one can only add spin (rotating or Kerr black holes) and charge (charged black holes) but no other 
independent characteristics: in the picturesque language used by Wheeler, it is said that black holes can 
have no hair. In a sense, the black holes of general relativity are very similar to fundamental particles, 
which are caracterized by a finite set of numbers (including mass, spin, and electric charge). 

Astrophysical black holes are somewhat more complex because of their material environment, as 
we will now see. 


3.3 Astrophysical black holes 


For a long time, black holes were considered as a curiosity of general relativity and did not have the status 
of other stellar objects. This has changed in the last decade which has seen mounting evidence that, at the 
center of our own galaxy (Milky Way) cluster, there is a massive black hole associated with the compact 
radio source Sagittarius A*. Observations of the motions of nearby stars by the imager/spectrometer 
NAOS/CONICA working in the infrared [46] have indeed confirmed the presence of a very massive 
object ((2.6 ± 0.2). 10 6 solar mass) localized in a very small region (a fraction of an astronomical unit, 
see Fig.[T2|), which seems only compatible with a black hole. 

Since then, black holes have been identified in many instances and their role might be central 
in many phenomena. We have seen that they are very simple gravitational objects. But these simple 
objects accrete matter, and are thus associated with very diverse phenomena. A picture has emerged, 
which seems to be valid at very diverse scales (see Fig. 13 1 of a black hole surrounded by an accretion 
disk and a torus of dust, with two opposite relativistic jets, which are supposedly formed during the 
gravitational collapse through a recombination of the magnetic fields. Of course, at the centre of this 
complex structure, lies the black hole surrounded by its horizon. But the complex phenomena that take 
place in this surrounding region allow to detect indirectly the black hole. 

Let us review some of the astrophysical occurences of black holes. 


First, our galaxy is not the only onewhich has a cventral black hole. This is believed to be very 
general, and in many cases the blck hole and its environment is much more active than our own. Active 
galaxies are galaxies where the dominant energy output is not due to stars. In the case of Active Galactic 
Nuclei (AGN), the non-thermal radiation comes from a central region of a few parsecs around the centre 
of the galaxy. The most famous example of such AGNs is provided by quasi-stellar objects (QSOs) or 
quasars: these starlike objects turn out to be associated with the point-like optical emission from the 
nucleus of an active galaxy. 
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Fig. 12: Left: NAOS/CONICA image of the central 2” (Sgr A* in light blue). Right: orbit of the star S2 around 
Sgr A* (4gj. 


The typology of active extragalactic objects is very complex: radio loud and radio quiet quasars, 
Seyfert galaxies, BL Lacs or blazars... There has been an effort to build a unified picture 1471: the 
apparent diversity in the observations would then result from the diversity of perspectives from which 
we observers see these highly non-isotropic objects. Typically, the model for radio-loud AGNs includes 
(see Fig. [13} right panel): a central engine, a pair of oppositely directed relativistic jets (cones of semi¬ 
angle around 1°) an accretion disk (of size of the order of 1 parsec), and a torus of material (of size of 
the order of 100 parsec) which obscures the central engine when one observes it sideways. Depending 
on the relative angle between the line of sight and the jet axis, observation may vary in important ways. 


Gamma ray bursts (GRB) are the most luminous events observed in the universe. They were 
discovered accidentally by the American military satellites VEGA which were designed to monitor the 
nuclear test ban treaty of 1963. The first burst was found in 1969, buried in gamma-ray data from 1967: 
two Vela satellites had detected more or less identical signals, showing the source to be roughly the same 
distance from each satellite £48]. 

A GRB explosion can be as luminous as objects which are in our vicinity, such as the Crab nebula, 
although they are very distant. The initial flash is short (from a few seconds to a few hundred for a long 
GRB, a fraction of a second for a short one). From 1991 to 2000, BATSE (Burst and Transient Source 
Experiment) has allowed to detect some 2700 bursts and showed that their distribution is isotropic, a 
good argument in favor of their cosmological origin. In 1997 (February 28), the precise determination 
of the position of a GRB (hence named 970228) by the Beppo-SAX satellite allowed ground telescopes 
to discover a rapidly decreasing optical counterpart, called afterglow. Typically in the afterglow, the 
photon energy decreases with time as a power law (from X ray to optical, IR and radio) as well as the 
flux: it stops after a few days or weeks. The study of afterglows gives precious information on the 
dynamics of GRBs. The launch of the SWIFT satellite on 20 November 2004 has started a new era for 
the understanding of GRBs. 

Given the time scales involved (a few milliseconds for the rise time of the gamma signal), the size 
of the source must be very small: it cannot exceed the distance that radiation can travel in the same time 
interval, i.e. at most a few hundred kilometers. Energy must have been ejected in an ultra-relativistic 
flow which converted its kinetic energy into radiation away from the source: the Lorentz factors involved 
are typically of the order of 100! 

This flow is collimated and forms a jet of half opening angle 6 . The observational evidence for 
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Fig. 13: Unified picture of the system surrounding an astrophysical black hole in the case of a gamma ray burst 
(left), a microblazar (centre) and a blazar (right). In the first two cases, the central object is a stellar mass black 
hole and the size of the accretion disk is about 1000 km. In the latter (AGN), it is a supermassive black hole of 
potentially several million solar masses and the accretion disk is of the order of 1 billion kms. 


this collimation is an achromatic break in the afterglow light curve: for t > f jeti , it decreases faster than 


it would in the spherical case [49 50]. If we assume that the relativistic jet, after emitting a fraction v/ 7 of 
its kinetic energy into prompt 7 rays, hits a homogeneous medium with a constant number density n, the 
break appears in the afterglow light curve when the Lorentz factor 7 becomes of the order of 1/6. This 


gives a relation between the half opening angle 6 and the break time fj et [50]: 


e = 


0.161 


(1 + z) 3 / 8 \1 day 


fjet 


3/8 


10 oz ergs 
E, 


7,iso 


1/8 


n 

1 cm 3 


1/8 


< /s 


(151) 


where i so is the isotropic equivalent gamma ray energy. 


In the collapsar model of Woosley [51 ], long GRBs are associated with the explosion of a rapidly 
rotating massive star which collapses into a spinning black hole (see Fig. [13] left panel). The burst and 
its afterglow have been successfully explained by the interaction of a highly relativistic jet with itself 
(internal shocks [52] 531) and with the circumstellar medium (external shocks |54|). Typically, one 
expects per day 10 6 collapses of massive stars in the Universe; 10 3 give rise to a GRB and approximately 
1 of these is pointing towards us its jet. Flence, one may observe from earth about one GRB per day. 


Supernovae explosions also provide very bright events in the sky, some of them being visible to 
the naked eye. Supernovae explosions were thus recorded in 1006, 1054, 1181, 1572 and 1604 . The Sn 
1987 A explosion allowed the detection of neutrinos and gamma emission. 
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The modem theory of supernovae was initiated in the 30s by Baade and Zwicky [551. 


Supernovae follow a classification according to spectroscopy. In type I supernovae, hydrogen 
lines are absent whereas they are present in type II. Moreover, type I has subclasses: for example, type 
la involves intermediate mass elements (Si). Each type corresponds to a different mechanism for the 
explosion. In particular, type II and type la have a completely different interpretation. 


We will focus first on type II supernovae. 


Presupernova stars (M > 8AT.,) have an onion-like structure. From the outer to the inner layers, 
one finds increasingly heavy elements: H, He, C, O, Ne, Si and Fe. 

As Si is consumed by nuclear reactions, the mass of the Fe core increases . The resulting density 
increase then turns the electrons relativistic and makes electronic capture (p + e — > n + v) energetically 
favorable. This diminishes the degenerate electron pressure and leads to the collapse of the core. Since 
Pcore ~ 10 12 kg.m -3 , the collapse time is typically ~ 0.1 s. 

This time, neutrinos produced as electrons are turned into neutrons, are trapped in the imploding 
core. The critical density for which neutrinos are trapped, is typically p ~ 2 x 10 14 kg/m 3 . As the core 
is crushed to higher densities, the density approaches that of a neutron star (p ~ 2 x 10 17 kg/m 3 ) and 
matter becomes almost incompressible. If the process was elastic, the kinetic energy would be enough 
to bring it back to the initial state. Typically 


E ~ G N M core 






G n M core 
Rns 


~ 3 x 10 46 J . 


(152) 


This is not completely so but there is a rebound of the core which sends a shock wave outward. 
Meanwhile, the stellar matter has started to free fall since it is no longer sustained by its core. The falling 
matter meets the outgoing shock wave and turns it into an accretion wave. 


Neutrinos emitted from the core heat up and expand the bubble thus formed. Convection and 
neutrino heating thus convey a fraction of the order of one percent of the neutron star gravitational mass 
(152 1 to the accretion front. This is enough to make it explode. 


One word of caution however: numerical models that try to reproduce supernovae explosions have 
been until now unable to explode the supernovae! One needs to start the explosion artificially. It therefore 
remains possible that one is still missing a key ingredient in the recipe. 


The bulk of the star blown off by the explosion makes what is known as a supernova remnant. 
It sweeps the interstellar medium at great velocity (10000 km/s) and may remain visible for 10 5 to 10 6 
years. A large fraction of the interstellar medium is thus swept by supernovae remnants (see Exercise 
3-1). This is important since this is believed to be the way the heaviest nuclear elements are scattered in 
the universe (primordial nucleosynthesis produces no element heavier than ' Li). 


Exercise 3-4 : a) Assuming approximately one supernova explosion every 30 years in our galaxy 
(assimilated to a disk of radius 15 kpc and thickness 200 pc), compute the corresponding rate 1Z of 
supernovae explosions per pc 3 and per year. 

b) If every supernova leads to a remnant of radius R = 100 pc that lasts for t ~ 10 6 yrs, what 
fraction of the galaxy volume is filled by the supernova remnant? 

Hints: a) TZ ~ 2.3 x 10 -13 pc~ 3 .yr _1 . 

b)l — exp [—(4-7r/3)ii 3 72i] ~ 0.5. 

SNIa events on the other hand are thermonuclear explosions of white dwarfs. More precisely, a 
carbon-oxygen white dwarf accretes matter (from a companion star or by coalescence with another white 
dwarf) which causes its mass to exceed the Chandrasekhar limit. The central core collapses, making the 
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carbon bum and causing a wave of combustion to propagate through the star, disrupting it completely. 
The total production of energy is thus almost constant. For a white dwarf of radius 1500 to 2000 km, 
about 2 x 10 51 ergs is released in a few seconds during which takes place the acceleration of the material. 
This is followed by a period of free expansion. Virtually all the energy of the explosion goes into the 
expansion. The luminosity of the supernova, on the other hand, finds its origin in the nuclear decay of 
the 56 Ni freshly synthesized. The energy release in the nuclear decays 56 Ni —> 56 Co —> 56 Fe, with 
respective lifetimes of 8.8 and 111 days, represents a few percent of the initial energy release. 


This model allows to understand the homogeneity of the observed type la supernovae explosions 
and why they have been used successfully as standard candles in cosmology (see Section 5.1 of Chap¬ 
ter [5]). The structure of a white dwarf is determined by degenerate electrons and thus independent of 
detailed chemical composition (see Section [XT] ). The rate of expansion is set by the total energy avail¬ 
able since the complete white dwarf is disrupted. Finally, the absolute brightness is determined by the 
radioactive decay of 56 M produced during the explosion. Less Ni means a lower luminosity but also 
lower temperature in the gas and thus lower opacity and more rapid energy escape. Thus dimmer super- 
novae are quicker i.e. have narrower light curves. 


3.4 High energy cosmic particles 

As explained above, compact objects and the violent phenomena associated with their formation are 
important to understand the origin of high energy cosmic particles. It is important to identify the potential 
sites of acceleration. Obviously, the jets described in the preceding Section are sources of energetic 
particles. Shock fronts, such as supernovae remnants, are also the siege of acceleration for particles 
whose multiple scattering off magnetized clouds lead to multiple encounters with the shock front. 

M. Hillas | [56| has proposed a general discussion of potential acceleration sites, in terms of the 
magnetic fields B availalble and the size R of the site. The Larmor radius of the particle rj j = E/(qBc) 
(in relativistic regime) may, with increasing energy E, become larger than the dimension R of the accel¬ 
erating site. We thus have the condition (q = Ze) 

E< Em = qB cR - Z (JL) (^) 9.3 10*>.V . (153) 

Note that this is the work of the electric field £ = Be over the maximal distance R. In the case of 
acceleration on magnetic clouds or a shock wave (where £ = BV), the maximal energy reads: 

E mSLX = ZeBVR . (154) 


In the case where acceleration involves large Lorentz factors 7, an extra factor 7 should be included to 
account for the energy being measured in the lab frame (also £ = 7 BV). 


This general criterion allows to draw the now classical Hillas diagram which identifies the possible 
acceleration sites in a plot log (B /1 G) vs. log (R /1 km). As can be seen on Fig. 14 given species of 
cosmic particles accelerated at given energies are represented by diagonal lines (from top to bottom on 
the figure: protons of 10 21 eV, protons of 10 20 eV and iron nuclei of 10 20 eV). 


4 Light does not say it all (2): dark matter 

We have known for a long time that the Universe has a dark component, and is not only luminous matter 
and radiation: already in 1933, by studying the velocity distribution of galaxies, Zwicky | [57| discovered 
that the Coma cluster had 400 times more mass than expected from its luminosity. Through the XX th 
century, it was realized that non-luminous matter, dark matter, is needed at all scales from galactic to 
cosmological. This pleads for a form of matter which is not included in the Standard Model, hence 
for physics beyond the Standard Model. Indeed, the detection of dark matter might be the first sign of 
physics beyond the Standard Model. Hence the programs of direct or indirect dark matter detection are 
of key importance not only for astrophysics but also for high energy physics. 
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Fig. 14: Hillas diagram showing size and magnetic fields of potential acceleration sites. Sites below the diagonal 
lines cannot accelerate protons above 10 21 eV, protons above 10 21 eV and Fe nuclei above 10 20 eV, respectively 
from top to bottom. 


4.1 The observational case 


As we alluded to above, dark matter was first identified by Fritz Zwicky |p7 58] in 1933 when studying 
the velocity distribution of galaxies in the Coma cluster. Using the virial theorem 


2(Uki n ) — —(E pot ) 


(155) 


(where (• • •) indicates time averaging), he concluded that there is 400 times more mass than expected 
from the luminosity. 

This was consistently confirmed by the study of the rotation curves of galaxies i.e. the velocity 
v(r) of stars as function of their distance r to the centre of the galaxy. Using again the virial theorem 
(155 l, we have for stars in the outer regions of the galaxy, since E y ~ rnv' 2 and E pot ~ G N mM/r (rn 
is the mass of the star, M of the galaxy), 


V oc 


G n M 


(156) 


Thus, one should see the velocity decrease as r 1//2 for stars at the border of the galaxyp^j This is not 
what is observed. Indeed, the rotation curves of galaxies, which were studied thoroughly through the 60s 


17 Note that for stars in the interior of the galaxy one should replace M by the mass of the sphere of radius r (M (r) oc p ss ,\r 3 ) 
which makes v(r) increase with r. 
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and 70s showed (see Figure [l5|) that the velocities do not start decreasing at the border of the luminous 
galaxy, as if there was more matter beyond. In fact, one needs a factor of order 10 more matter in spiral 
galaxies. The extra matter forms a halo that extends beyond the luminous galaxy. 


DISTRIBUTION OF DARK MATTER IN NGC 3198 



Fig. 15: Rotation curve v{r) for the galaxy NGC 3198. The curve labeled “disk” indicates the curve due to the 
stars in the galaxy, which extend only to 10 kpc; the curve labeled “halo” is the one that would be due solely to a 
spherical halo of dark matter. 


For example, the modem picture of our own Galaxy, the Milky Way, is one of a luminous bulge 
of a few kpc at the centre of a disk of radius 12.5 kpc and thickness 0.3 kpc, containing some 10 11 stars, 
surrounded by a nearly spherical halo of dark matter of typical radius 30 kpc (see Appendix [A] for the 
definition of a parsec). 

But dark matter is not only present in galaxies. We have seen that it was first identified by Zwicky 
in clusters of galaxy. This is now confirmed by many observation of clusters. X-ray studies have revealed 
the presence of large amounts of intergalactic gas which is very hot, and hence emits X-rays. The total 
mass of the gas is greater than that of the galaxies by roughly a factor of two. However this is not enough 
mass to keep the galaxies within the cluster. Since this gas is in approximate hydrostatic equilibrium 
with the cluster gravitational field, the mass distribution can be determined, which leads to a total mass 
estimate approximately six times larger than the mass of the individual galaxies or of the hot gas. 


A powerful tool for mapping dark matter is gravitational lensing, which is based on the deflection 
of light by matter: the light of a distant galaxy is deflected by an accumulation of matter present on the 
line of sight, just as it is by a lens (see Fig. 16). The deviation of light rays depends on the ratio of 
distances between observer, lens and source, as well as on the mass of the deflector. 

More precisely, the “lens equation” may be with written as (see Fig. [l5]for notations): 


0i 


a _i_ Dls 

Os + jx —« 
T'OS 


(157) 


It is well-known that the deflection angle of a light ray passing an object of mass M with an impact 
parameter b is a = 4 G N M/{bc 2 ) {see for example / 59 / p.286). Hence the lensing effect allows to detect 
the mass distribution. 
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Fig. 16: Configuration of a gravitational lens: S is the astrophysical source, O is the observer, L is the (point) lens 
and I is the image. The source lies at an angle 9g from the axis OL but is seen at an angle 9j: the deviation angle 
is thus a. 


Massive clusters may induce multiple images of background galaxies. One then talks of strong 
lensing [60|. The shape of a single galaxy can also be deformed into an arclet through lensing. In the case 
of weak lensing, the effect is measured through the deformation of the shape of galaxies but, because 
galaxies do not have a circular shape, it can only be measured statistically: galaxies tend through lensing 
to have aligned shapes j5TJ. 

Finally, we have seen in Section [13] that cosmological data shows that dark matter is also needed 
at the largest scales. An illustration of this is the map of dark matter that the Planck collaboration could 
draw using the gravitational lensing of the CMB light (see Fig.[T7]). 



Fig. 17: Large scale structures detected by Planck through their lensing potential (galactic North) 


Since the only proof of existence of dark matter is gravitational (rotation curves, lensing,...), one 
may wonder whether the observed phenomena are due to a modification of gravity, which would then 
be different from what general relativity predicts at the corresponding scales. Besides the difficulty of 
finding a theory that encompasses all the successes of general relativity, the problem is to modify gravity 
at all the scales where we see signs of dark matter, that is galaxies, clusters of galaxies and cosmological 
scales.For example, the MOND theory [63j64| has been proposed to explain the rotation curves of 
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galaxies but is only Newtonian and requires to be generalized [65] in order to be valid at the scale of the 
Universe. 


The existence of the bullet cluster (see Fig. 181 where two galaxies collide has been presented as a 
support of dark matter 1661 because the luminous parts of the galaxies are displaced with respect to their 
halos: because dark matter is weakly coupled, the halos (detected through gravitational lensing) continue 
their way during the collision whereas their - luminous matter counterparts (detected through their X-ray 
emission) are deformed. 


CD 

LO 


m 


oo 

m 

°tn 

LO 

I 


6 h 58 m 42 s 36 s 30 s 24 s 18 s 12 s 


Fig. 18: X-ray image of the merging cluster 1E0657-558 obtained by the satellite Chandra, over which is super¬ 
imposed in green contours the weak lensing reconstruction of the dark matter halos | 66) . 



4.2 Dark matter particles 

What is dark matter? Since it is nonluminous it is not electrically charged, and the only possible candidate 
within the Standard Model is the neutrino. But the random motion of neutrinos (the technical term is “free 
streaming”) would wash out any density fluctuation and prevent the formation of galaxies; one expresses 
this by saying that neutrinos are hot dark matter. Instead, we need cold dark matter i.e. particles with 
smaller free streaming length. 

Moreover, we need dark matter particles in sufficient quantity. This means that they cannot be 
in thermal equilibrium today: they must have decoupled from the thermal history of the Universe at 
some early time. Typically, there are two competing effects to modify the abundance of a species X: 
XX annihilation and expansion of the Universe. Indeed, the faster is the dilution associated with the 
expansion, the least effective is the annihilation because the particles recede from one another. When 
the temperature drops below the mass m x . the annihilation rate becomes smaller than the expansion rate 
and there is a freezing of the number of particles in a covolume. In more quantitative term, this reads for 
the freezing temperature Tf. 

n X ( T f) < a annV >~ H{Tf) , (158) 


where < er ann n > is the thermal average of the XX annihilation cross-section times the relative velocity 
of the two particles annihilating. One finds for the present density (in units of p c as usual) (see Ref. /[67j/ 
Section 5.5) 


n x K 


2 1.07 x 10 9 GeV’ 1 

rs j - 


Xf 


(159) 


where Xf = m x /(kTf ~ 20 and p* is the total number of relativistic degrees of freedom present in the 
universe at the time of decoupling. 
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We note that the smaller the annihilation cross section is, the larger is the relic density. We find 
n.Y ~ (100 TeV) -2 ( < a a n„ v > ) _1 ~ 0.1 pb/ < a ann v > (in units where f, = c = 1, 1 pb 
= 2.5 x 10 -9 GeV~ 2 ). This should be compared with the latest result Qdm = 0.1187 ± 0.0017 coming 
from Planck [16}. 

Thus Q x will be of the right order of magnitude if < a :inn v > is of the order of a picobarn, which 
is a typical order of magnitude for an electroweak process. Also writing dimensionally 

a 2 

< cr ann v > ~ — 2 ~ , (160) 

m x 


where a is a generic coupling strength, we find that is of order 1 for a mass mx ~ o x 1000 TeV, 
i.e. in the TeV range. This is why one is searching for a weakly interacting massive particle (or wimp). 

There is also the possibility that dark matter particles have produced non-thermally, e.g. from the 
decay of heavy particles. 


A puzzle which is not addressed by the wimp scenario is why dark matter and baryonic matter 
densities are basically of the same order: 

DrM\/T 

~ 5 . (161) 


Pdm 


PB 

Indeed, baryon density arises from baryogenesis (see section [L6| > and thus results from a small mismatch 
between baryons and antibaryons, as seen from (58). On the other hand, the wimp density results from 
the freezing regime described by (158). There is no reason that the two scenarios lead to similar energy 
densities, as in (161). This puzzle is addressed by the Asymmetric Dark Matter scenarios (see the review 
by C. Zurek [68] and references therein). The idea is that dark matter has an asymmetry in the number 
density of matter over anti-mattei[^] similar to the one for baryons: 


n x -nx~n b -n- b . 


(162) 


The abundance is therefore approximately one part in 10 iU in comparison with the thermal abundance 
(see (581). Eq. (162) suggests that mx is typically 5 times the proton mass, as a typical baryon mass. 
Thus generic Asymmetric Dark Matter models tend to favor light dark matter particles. 


One may search for dark matter particles through direct detection using their elastic collisions with 
nuclei XN —> XN in ultra-low background detectors. The energy of the recoiling nucleus is typically 
from a few keV to tens of keV. The recoil rate after integration over the dark matter velocity v distribution 
is 


R ~ 


3.5 x 1(T 2 


events 100 
kg.day A 


lOOGeV &xn (v) p 

m x X lpb X 220km.s _1 X 0.3GeV.cm~ 3 J ’ 


(163) 


where A is the atomic mass of the recoil nucleus, <jxn the cross-section for dark matter particle-nucleus 
elastic scattering and p the local density of dark matter in our Galaxy. In the case of a neutralino wimp, 
(Txn can be as low as 10 12 pb, which yields a rate of 10 -8 events/ton.year! Present, and future, 
experimental limits are shown on Fig. 19 One should note that there is below 10“ 12 pb (10“ 8 pb 
for low-mass particles) an irreducible neutrino background corresponding to the reaction vN —> uN. 
Experiments of the next decade should be able to reach this limit. 


Dark matter particles may also be searched for through their annihilation products in massive 
celestal objects. This is known as indirect detection. Indeed, because they are massive, they tend to 
accumulate in gravitational potentials, such as the centre of the Sun, or the centre of our Galaxy. They 
annihilate there into pairs of energetic particles, the energy of which is directly connected with the mass 
mx of the dark matter particles. One may therefore search for excess of energetic particles (positrons, 


18 This obviously precludes models where the dark matter particle is its own antiparticle, as often the case for wimps (e.g. the 
neutralino). 
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Fig. 19: Sensitivities of some running and planned direct detection dark matter experiments to the spin-independent 
elastic scattering cross-section. Full curves correspond to limits from existing experiments, dashed curves to 
predicted sensitivities of future experiments. The full brown, pink, blue and yellow regions correspond respectively 
to the regions allowed by potential signals observed by the DAMA, CREST, CDMS and CoGeNT experiments. 
The light red disk corresponds to the region favoured by supersymmetric models. The thick yellow line corresponds 
to the irreducible neutrino background. 


photons, neutrinos) in the direction of galactic centres such as in our own Milky Way. An excess of ener¬ 
getic positrons (energy of a few tens to few hundred GeV) has actually been observed by the PAMELA, 
Fermi and AMS-02 experiments. It remains to be seen if this arises from the annihilation of dark matter 
or from astrophysical sources, such as pulsars: we have seen at the end of Section[3]that a certain number 
of astrophysical sources produce energetic particles. This is indeed a limitation, at least at present, of the 
indirect detection of dark matter: it needs to be complemented by either direct detection or detection at 
colliders. 

4.3 WIMPs and physics beyond the Standard Model 

We would like to stress in this section that the presence of a WIMP in a theory is deeply connected with 
the naturalness of the electroweak scale. 

Let us start by recalling what is the naturalness problem (see for example l~67]j). As is well-known, 
the Higgs squared mass m 2 receives quadratically divergent corrections. In the context of an effective 
theory valid up to a cut-off scale A where a more fundamental theory takes over, A is the mass of the 
heavy degrees of freedom of the fundamental theory. Their contribution in loops, quadratic in their 
mass, destabilizes the Higgs mass and thus the electroweak scale (ml ~ Xv 2 where A is the scalar self¬ 
coupling and v i /(G F V2y/ 2 250 GeV is the Higgs vacuum expectation value. More precisely, we 


44 







have at one loop 


9 3 ml 

dmt = 


A 2 - 
2tt' 2 v 2 H 


6MI+3M 2 

8tt 2 v 2 9 


3 m l 

8 n 2 v 2 


A 2 


h > 


(164) 


where for completeness we have assumed different cut-offs for the top loops (At), the gauge loops (A g ) 
and the scalar loops (A/J The naturalness condition states that the order of magnitude of the 

Higgs mass is not destabilized by the radiative corrections i.e. dm 2 . < m 2 . This translates into the 
conditions: 


A t ~ 

Ag ~ 

Ah ~ 


\ — m h ~ 8.5m h , 
3 rn t 

2 \f2irv 


f M l + Z M l 

2 \J~2 , kv 


m h ~ 9m h , 


V3 


1.3 TeV. 


(165) 

(166) 

(167) 


Thus one should introduce new physics at a scale At ~ 3.5?/)/,. We will illustrate our argument with 
two examples: supersymmetry and extra dimensions. In the two cases, one introduces new physics at the 
scale At (supersymmetric particles or Kaluza-Klein modes). 

Typically, these models require the presence of a symmetry that prevents direct coupling between 
the Standard Model (SM) fermions and the new fields that one has introduced: otherwise, such couplings 
introduce new mixing patterns incompatible with what is observed in flavor mixings (compatible with the 
Standard Model). This symmetry is usually a parity (i.e. a discrete symmetry) which is the low energy 
remnant of a continuous symmetry which operates at the level of the underlying fundamental theory: SM 
fermions are even under this parity whereas the new fields are odd. Among these new fields, the lightest 
odd-parity particle (we will refer to it as the LOP) is stable: it cannot decay into SM fermions because 
of the parity; it cannot decay into the new fields because it is the lightest. It is massive and weakly 
interacting. It thus provides an adequate candidate for a WIMP. 

Let us take our examples in turn. In the case of supersymmetry, the parity operation is R-parity 
(which usually proceeds from a continuous R-symmetry broken by gaugino masses i.e. supersymmetry 
breaking). And the LOP is the Lightest Supersymmetric Particle, the famous LSP, the lightest neutralino 
in the simplest models. 


In the case of extra dimensions, say a 5-dimensional model, the local symmetry is 5-dimensional 
Lorentz invariance. It ensures conservation of the Kaluza-Klein levels: if is the nth Kaluza-Klein 
mode of the massless 5-dimensional field A (in other words, the ^-dimensional field with mass m = n/II, 
where R is the radius of the 5th dimension), then in the reaction A + B^ p ) —> + D^ r \ we have 

n + p = q + r. At energies smaller than R~ l , this turns into a Kaluza-Klein parity (—l) n . The LOP is 
then the lightest Kaluza-Klein mode, usually If 1 the first mode of the U{\)y gauge boson / 70 /. 


In realistic models, there is often the possibility that other odd-parity fields are almost degenerate 
in mass with the LOP. This leads to the possibility of co-annihilations, that is annihilations of the LOP 
against these almost degenerate fields, and to a modification of the relic density in the corresponding 
region of parameter space. 


Searches at LHC are based on the missing energy signal corresponding to the LSP (see Fig. 20 1 . 
Since LSP are produced in pairs, they are difficult to reconstruct in all generality. But, in the case of a 
specific model, one may be able to reconstruct the mass of the LSP as well as the relic density. 


4.4 Other candidates for dark matter 

Among the many other candidates proposed for dark matter, one may single out the axion field since it is 
introduced to solve one puzzle of the Standard Model known as the strong CP problem, which remains to 
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Fig. 20: Upper panel: Contours in a parameter space of supersymmetry models in the plot (mo, m\/ 2 ) for the 
discovery of the missing energy plus jets signature of new physics by the ATLAS experiment at the LHC. The 
three sets of contours correspond to levels of integrated luminosity at the LHC (in fb -1 ), contours of constant 
squark mass, and contours of constant gluino mass (ZD- 

Lower panel: Exclusion limits at 95% CL for 8 TeV, 20 fb _1 integrated luminosity analyses in the (mo, rrii/ 2 ) 
plane for the MSUGRA/CMSSM model with the remaining parameters set to tan/? = 30, A 0 = —2?n 0 , H > 0. 
Part of the model plane accommodates a lightest neutral scalar Higgs boson mass of 125 GeV fZD- 


be solved. Moreover, from the point of view of cosmology, this is an interesting illustration of a low-mass 
and wery weakly interacting particle, as we will encounter in the next Chapter. 

As we have already seen, CP symmetry is violated in weak interactions. But the symmetries of the 
Standard Model allow as well a CP-violating term in the QCD Lagrangian (#p 2 /327r 2 )G'“ l/ G aAt % where 
is the gluon field strength and g the QCD gauge coupling. For non-zero quark masses, this term 
leads to unobserved CP-violating effects in the strong sectoip*) 

The most common way to solve the puzzle is to introduce a scalar field a(x) called axion, with 
Lagrangian 

C a = \d»ad,a + , (168) 


where f a is an energy scale, called the axion decay constant (by analogy with the pion decay constant). 
Strong interactions generate an effective potential for a(x) whose minimum corresponds to no violation 
of Cf} 2 °! The non-renormalisable interaction aG ■ G is obtained as the low energy effect of the breaking 
of a (7(1) global symmetry, known as the Peccei-Quinn symmetry [73 74], spontaneously broken at the 
scale f a : the axion is the pseudo-Goldstone boson associated with this breaking [75 76 pM 


19 To be more precise, it is 9 = 9 — arg det m q , where m q is the quark mass matrix, which is obsen’able. If 9 A 0, then 
strong interactions violate P and CP. This is not compatible with the experimental upper bound on the neutron electric dipole 
moment unless |#| < 10 -10 . 

20 i.e. a vanishing value of 9 = a(x)/f a — arg det m q 

21 The axion is not a true Goldstone boson but a “pseudo-Goldstone" boson because the QCD vacuum which involves non¬ 
vanishing values for the quark condensates such as (uu) and {dd ) indices a small explicit breaking of the U( 1) symmetry, 
hence a mass for the otherwise massless Goldstone boson associated with the spontaneous breaking. 
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The axion mass is 


10 12 GeV 


m a ~ 6.10~ 6 eV 


fa 


(169) 


Its coupling to ordinary matter is proportional to 1 / f a and can be calculated in specific models. It couples 
to leptons and to photons, the latter being of the form 


£ 


a 77 — 


9a-y 


a a(x) 

vr fa 


E.B , 


(170) 


where E and B are the electric and magnetic fields, a is the fine structure constant and g ai is a model- 
dependent coefficient of order 1. Moreover, since m a -C A qcd, the axion coupling to quarks should be 
described through its coupling to hadrons, which occurs through small mixing with the ttq and // mesons. 
All of these interactions can play a role in searches for the axion, and allows the axion to be produced 
or detected in the laboratory and emitted by the sun or other stars. Its non-discovery leaves us with an 
axion window 10 -6 eV < m a < 3.ICE 3 eV, or correspondingly, 2.10 9 GeV < f a < 6.10 12 GeV. 


Let us describe briefly the cosmology of the axion field (see the review by P. Sikivie [77] for a 
more thorough treatment). The breaking of the (7(1) symmetry corresponds to a phase transition, known 
as the Peccei-Quinn phase transition, at a temperature of order f a . This phase transition is characterized 
by the formation of comic strings. 


If the reheat temperature after inflation is smaller than f a , then one starts the evolution in the 
reheated universe with an homogeneous axion field. When the temperature reaches the QCD scale, the 
effective potential turns on and the axion acquires a mass. At a time i* ~ m a 1 ( T * 1 GeV), the axion 

starts to oscillate around its minimum. These oscillations do not dissipate into other forms of energy, and 
thus contribute to the cosmological energy density an amount 





0.15 


( fa 


V 10 12 GeV 


7/6 



(171) 


where a(t*) is the axion value at t *, which measures the misalignment of the axion with respect to its 
final minimum. Such a contribution is thus called vacuum realignment. 

If the reheat temperature is larger than f a , then the cosmic strings produced at the Peccei-Quinn 
phase transition become at time £ the boundaries of domain walls. This leads to a potential domain wall 
problem (too much energy stored in the domain walls). There is a certain number of cases where this can 
be avoided. In this case Gz)> 





0.7 


( fa 


\ 10 12 GeV 


7/6 


(172) 


5 Light does not say it all (3): dark energy 

We have seen in the introduction that the observation of the acceleration of the expansion of the Universe 
in 1998-1999 provided a way out of an increasingly uncomfortable tension between models and obser¬ 
vation. On one hand, the observed (luminous and dark) matter could account for only a fraction of the 
critical density p c = 10~ 26 kg/nr 3 (other forms of energy being negligible at present time). On the other 
hand, the standard theory of inflation erased any curvature and naturally led to a spatially flat universe 
for which p = p c . The discovery of the acceleration of the expansion led to introduce a new component, 
named dark energy, of a type unknown so far since all known forms of energy (non-relativistic matter, 
radiation) decelerate the expansion. 

In this last chapter, after briefly reviewing the observational case, we will illustrate dark energy 
models with the example of quintessence, identify some of the problems posed by dark energy models, 
and the fundamental questions associated with the acceleration of the expansion of the Universe. 
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5.1 Acceleration of the expansion of the Universe: supernovae of type la as standard candle 

The approach that has made the first case for the acceleration of the expansion of the Universe uses 


supernovae of type la as standard candles (see Section 3.31. Two groups, the Supernova Cosmology 


Project [7] and the High -2 Supernova Search [6j have found that distant supernovae appear to be fainter 
than expected in a flat matter-dominated Universe. If this is to have a cosmological origin, this means 
that, at fixed redshift, they are at larger distances than expected in such a context and thus that the 
Universe expansion is accelerating. 

More precisely, one uses the relation (C.31 between the flux d> received on Earth and the luminosity 


L of the supernova. Traditionally, flux and luminosity are expressed on a log scale as apparent magnitude 
ms and absolute magnitude M (magnitude is —2.5 log 10 luminosity + constant). The relation then reads 


= 5 log (HocIl) + M — 5 log H 0 + 25. 


(173) 


The last terms are ^-independent, if one assumes that supernovae of type la are standard candles; they 
are then measured by using low 2 supernovae. The first term, which involves the luminosity distance 
varies logarithmically with 2 up to corrections which depend on the geometry, more precisely on go for 


small 2 as can be seen from (C.6). This allows to compare with data cosmological models with different 


components participating to the energy budget, as can be seen from Fig. 21 



Fig. 21: Hubble plot (magnitude versus redshift) for Type la supernovae observed at low redshift by the Calan- 
Tololo Supernova Survey and at moderate redshift by the Supernova Cosmology Project. 


U, 


In the case of a model with matter and cosmological constant as dominant components, 
f /2 — U,\ and the measurement can be turned into a limit in the Qm — -h\ plane see Fig. 22). 


Fet us note that this combination Q M /2—Q\ is ‘orthogonal’ to the combination 1—= Qm+^a 
measured in CMB experiments. The two measurements are therefore complementary: this is sometimes 
referred to as ‘cosmic complementarity’. 


An important question raised by the analysis above is whether supernovae are truly standard can¬ 
dles. Otherwise, the observation could be interpreted as an history effect: for some reasons, older su- 
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Supernova Cosmology Project 
Perlmutter etal. (1998) 



Fig. 22: Best-fit coincidence regions in the f 1m — U\ plane, based on the analysis of 42 type la supernovae 
discovered by the Supernova Cosmology Project 0 - 


pernovae would be dimmer. Indeed, strictly speaking, supernovae of type la are not standard candles: 
dimmer supernovae are quicker (see Section |d. J[ ). In practice, one thus has to correct the light curves 
using a phenomenological stretch factor. It is thus more precise to state that supernovae of type la are 
standardizable candles. Moreover, the type of measurement discussed above is sensitive to many pos¬ 
sible systematic effects (evolution besides the light-curve timescale correction, presence of dust, etc.), 
and this has fuelled a healthy debate on the significance of supernova data as well as a thorough study of 
possible systematic effects by the observational groups concerned. 


5.2 Cosmological constant and vacuum energy 

An obvious solution to the acceleration of the expansion is the introduction of a cosmological constant 
(see (451 remembering that go is a deceleration parameter). But there is a severe conceptual problem 
associated with the cosmological constant, which we now describe. 

Considering ([25]) in flat space at present time implies the general following constraint on A: 


|A| < Hi . 


(174) 


In other words, the length scale = |A| 1,/2 associated with the cosmological constant must be larger 
than the Hubble length Ih 0 = cH^ 1 = /iq 1 .10 26 m, and thus be a cosmological distance. 

This is not a problem as long as one remains classical: £h 0 provides a natural cosmological scale 
for our present Universe. The problem arises when one tries to combine gravity with the quantum theory. 
Indeed, from Newton’s constant and the Planck constant h, we can construct the (reduced) Planck mass 
scale 

m p = \Jhe/(8 itG n ) = 2.4 x 10 18 GeV/c 2 . (175) 


49 








The corresponding length scale is the Planck length 


l 


p 


h 

m p c 


= 8.1 x 10 _35 m . 


(176) 


The above constraint now reads: 


i\ = iai - 1 / 2 >i Ha = 4~- 1060 4 • 

-Wo 


(177) 


In other words, there are more than sixty orders of magnitude between the scale associated with the 
cosmological constant and the scale of quantum gravity. 

A rather obvious solution is to take A = 0. This is as valid a choice as any other in a pure gravity 
theory. Unfortunately, it is an unnatural one when one introduces any kind of matter. Indeed, set A to 
zero but assume that there is a nonvanishing vacuum (i.e. ground state) energy: < T jnj >= p V ac <7/^: then 
the Einstein equations ([3]) read 

Ftfiv — 8ttG n T^ u -T 8ttG n p-va,cSfu/ ■ (178) 


As first noted by ZeTdovich [781, the last term is interpreted as an effective cosmological constant (from 
now on, we set h = c = 1): 

A 4 

A eff = 8 irG N p va , c = —3- . ( 179 ) 

m p 

Generically, p v;ic receives a non-zero contribution from symmetry breaking: for instance, the scale A 
would be typically of the order of 100 GeV in the case of the electroweak gauge symmetry breaking 
or 1 TeV in the case of supersymmetry breaking. Moreover, it is divergent in the context of the (non- 
renormalizable) theory of gravity, which would thus favour a value as large as the Planck scale. But the 
constraint ( |177[ ) now reads: 

( 180 ) 


A < 10 -30 m p ~ 


10" 3 eV. 


It is this very unnatural fine-tuning of parameters (in explicit cases p mc and thus A are functions of the 
parameters of the theory) that is referred to as the cosmological constant problem, or more accurately the 
vacuum energy problem. 

If the acceleration observed is indeed due to the cosmological constant, its value is as large as the 
upper bounds obtained in the previous subsection allow: 


A ~ Hq , ~ Ih 0 , A ~ 10 3 eV . (181) 

Regarding the latter scale A, which characterizes the vacuum energy (p v ac = A 4 ), one may note the 
interesting numerical coincidence: 


J ~ Wp ~ 10“ 4 m . ( 182 ) 

This relation underlines the fact that the vacuum energy problem involves some deep connection between 
the infrared regime (the infrared cut-off being £h 0 ) and the ultraviolet regime (the ultraviolet cut-off being 
£ p ), between the infinitely large and the infinitely small. 

As an illustration of what such a relation could tell us about some very fundamental aspects of 
physics (if it is not a mere numerical coicidence), we will follow the approach developed by T. Padman- 
abhan [79]. As we will discuss later, the value of the cosmological constant might be related to the 
way that the short distance theory reacts to long distance fluctuations. Let us consider a 3-dimensional 
domain of size L (e.g. the horizon £h 0 )- From the point of view of the quantum theory, it consists of 
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N = ( L/t p ) 3 elementary cells. For each individual cell, a “natural value” for the energy stored is 
provided by the scale m p characteristic of quantum gravity. This yields 


p ~ 



1 



(183) 


which is some 120 orders of magnitude larger than observed. 

Alternatively, if a mechanism, yet to be determined, cancels this bulk energy, vacuum energy may 
be produced by the energy fluctuations. The Poissonian fluctuation in energy is Ae ~ 1/^ P > which 
corresponds to an energy for overall fluctuations AE 2 ~ N/i 2 p , or an energy density 


Vn 


(184) 


" e p o fn L w ’ 

which again does not reproduce ( 182 ) fp vac = A 4 ). 

If we make the further assumption that the relevant degrees of freedom lie on the surface (as the 
degrees of freedom of a black hole lie on the horizon), then N = ( L/i p ) 2 and 


Vn 


i 


,L 3 i 2 L 2 ’ 


(185) 


which is fidly consistent with (182). Is this telling us something on the quantization of spacetime, hence 
of gravity? We will return to this question below. 


5.3 Supersymmetry 


The most natural reason why vacuum energy would be vanishing is a symmetry argument. It turns 
out that, among the various spacetime symmetries available, global supersymmetry is the symmetry 
intimately connected with the vanishing of the vacuum energy. 

We recall that supersymmetry is a symmetry between bosons and fermions which plays an im¬ 
portant role in high energy physics, mainly because, through a cancellation between the boson and the 
fermion fields, it controls severely the quantum fluctuations. Many believe that the Standard Model is 
the effective theory of a more fundamental theory valid at higher energies, with more boson and fermion 
fields. These fields should, to some level, make themselves known through the quantum energy fluc¬ 
tuations to which they participe. The fact that we find no trace of them seems to suggest that these 
fluctuations are tightly constrained, i.e. that the underlying theory is supersymmetric. 

Among the various spacetime symmetries, supersymmetry is rather unique. Indeed, in the same 
way as the generator of time and space translations is the 4-momentum operator P^ = (P° = //. P '), 
and the generator of spacetime rotations is the tensor M^ v , the only other generators carrying a Lorentz 
index are the generators of supersymmetry Q r , where r is a spinor index. In fact, the combination of two 
supersymmetry transformations is merely a translation in spacetime. This is expressed by the following 
algebra: 

{Qr: Qs} = H l s P n , 086 ) 

where Q = Qy° and -f 1 are the gamma matrices introduced by Dirac to write a relativistic fermion 
(electron) cquat iorp’j Since the generator of time translations Pq is the Hamiltonian //, we may easily 
infer from ( |l86[ ) an expression for the Hamiltonian of the systenj^j 

H = \Y^Ql • (190) 

r 


z The anticommutator in (186) arises from the fact that the supersymmetry transformation parameter is an anticommuting 
spinor. 


3 Indeed, {186) reads explicitly 


{Qr, Qt} 7i°s = ^rs Pli ■ 


( 187 ) 
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It follows that the energy of the vacuum |0) can be expressed as: 


<0|tf|0> = ^|| Q r \0) 


( 191 ) 


Thus, the vacuum energy vanishes if and only if supersymmetry is a symmetry of the vacuum: Q r \0) = 0 
for all rW\ 

The problem however is that, at the same time, supersymmetry predicts equal boson and fermion 
masses and therefore needs to be broken since this is not observed in Nature. The amount of breaking 
necessary to push the supersymmetric partners high enough not to have been observed yet, typically A ~ 
TeV, is incompatible with the limit ( 1 8Q[ ). 

Moreover, in the context of cosmology, we should consider supersymmetry in a gravity context 


and thus work with its local version, supergravity (following (1861, local supersymmetry transformations 
are associated with local spacetime translations which are nothing but the reparameterizations which 
play a central role in general relativity). In this context, the criterion of vanishing vacuum energy is 
traded for one of vanishing mass for the gravitino, the supersymmetric partner of the graviton (which 
allows to cancel the constant vacuum energy at the expense of generating a mass m 3/2 for the gravitino 
field; see e.g. Ref. [ 571 , section 6.3 for a more complete treatment). Local supersymmetry is then 


absolutely compatible with a nonvanishing vacuum energy, preferably a negative one (although possibly 
also a positive one). This is both a blessing and a problem: supersymmetry may be broken while the 
cosmological constant remains small, but we have lost our rationale for a vanishing, or very small, 
cosmological constant and fine-tuning raises again its ugly head. 

In some supergravity theories, however, one may recover the vanishing vacuum energy criterion. 


5.4 Why now? 

In the case where the acceleration of the expansion is explained by a cosmological constant, one has 
to explain why this constant contribution appears to start to dominate precisely now. This is the “Why 


now?” or cosmic coincidence problem summarized in Fig. 23 the coincidence being between the onset 
of acceleration and the present time (on the scale of the age of the Universe). In order to avoid any 
reference to us (and hence any anthropic interpretation, see below), we may rephrase the problem as 
follows. Why does the dark energy starts to dominate at a time t\ (redshift z\ ~ 1) which almost 
coincides with the epoch to (redshift zq ~ 3 to 5) of galaxy formation? 


Exercise 5-1 : In this exercise, we will study the evolution of a flat (k = 0) universe with a non- 
relativistic matter component (energy density p M ) and a cosmological term (or equivalently vacuum 


energy density p\ = X/(8 ttG n )). The relevant equations of evolution are obtained from (23), (29) and 


(31): 


d 2 


TT - 87T G n 

H - ^2 “ —3— (Pm + fa) > 


a 

a 


4ttG n , , 

2 \Pm 5 


Contracting with 7 sr , one obtains 


E Qt} [(7°) 2 ] tr = 2 TV ( 7 V) P 


(192) 

(193) 

(188) 


Using (7°)2 = 1 and Tr( 7 ° 7 M ) = 4<? , one obtains 


E Ql = 4 P° = 4 H . 


( 189 ) 


24 Remember that a supersymmetry transformation U is obtained by exponentiating the generators: t/|0) = |0). 
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Fig. 23: The cosmic coincidence problem illustrated in the case of a cosmological constant 


Pm = -3Hp M ■ 


(194) 


a) Defining H\ = (8itG n pa/3 ) l ,/2 , show that the solution to the system of equations is given by 


a{t) 


C (sinh -H\tj 


\ 2/3 


(195) 


where C is a constant. Compute the Hubble parameter H as well as the ratio p M /p\ as a function of 
time t and H\. 

b) Compute the present Hubble constant Ho in terms of H\ and Q v; . What is the age of the 
Universe to in terms of Ho and it M ? 

c) Plot a(fi)/C and P m /pa hi terms of time. When is the equality between matter and vacuum 
energy reached? 


Hints: a) H = H\/ tanh(3fVAf/2) and p M /p\ = [s'mh(fiH t \t/2)] 2 . 
b) H 0 = Ha/^T^T^, 

2 1 1 + y^^U 

0 3H 0 yr^; n ^ ■ 


(196) 


5.5 More dynamics: dark energy vs modification of gravity 

An alternate possibility is that the cosmological constant is much smaller or even vanishing and that 
the acceleration is due to some new form of energy -known as dark energy- or some modifications of 
gravity. These two possibilities correspond to modifications of either sides of Einstein’s equations ([3]). 
Let us envisage briefly these two cases. 

First, we may try to identify a new component p x of the energy density with negative pressure: 

Px= w xPx , W X < 0 • (19V) 
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Note that the equation of state parameter w x may not be constant and may thus evolve with time. 


Observational data constrains such a dynamical component, referred to in the literature as dark 
energy, just as it did with the cosmological constant. For example, in a spatially flat Universe with only 
matter and this unknown component X, one obtains from ( |3Tj ) with p = pm + p x , p = w x Px the 
following form for the Hubble parameter and the deceleration parameter (compare with (43 I and (|46j)) 


H\z ) 

Q(z) 


m 


\l M (l + zf + Sl x (l + zf^ w ^ 


Hi 


2 H(z) 


^m(1 + z ) 3 + + 3m x )(l + 2:) 3( - 1 +u 'x) 


where Ox = p x /p c . At present time (compare with (45 l), 

Qo — -pr- + (1 + 3w x ) — 


(198) 

(199) 


( 200 ) 


The acceleration of the expansion observed requires that Q x dominate^] with w x < —1/3. 


In Fig. 24 we present constraints in the (f^ 


plane, obtained recently [ 80 ] from combining 


observations using supernovae and BAO as well as CMB results from WMAP. 


-0.5 



Fig. 24: Confidence contours in the (f l M ,w x ) plane arising from supernova results from SuperNova Legacy 
Survey (SNLS) 3 year results (in blue) and combined BAOAVMAP7 constraints (in green) [see Ref. |80)]. 


An important property of dark energy is that it does not appear to be clustered (just as a cosmo¬ 
logical constant). Otherwise, its effects would have been detected locally, as is the case for dark matter. 
This points towards scalar fields, which generically have this property. Indeed, an attractive property of 
scalar fields is that they easily provide a diffuse background by resisting gravitational attraction. The key 

25 One may easily obtain from \l99\ the time of the onset of the acceleration phase: 


1 + Zacc — 


— (1 + 3w x ) 


-l/(3tu v ) 


( 201 ) 
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quantity when discussing gravitational clustering is the speed of sound defined as 




5p 

dp 


( 202 ) 


It is a measure of how the pressure of the field resists gravitation clustering. In most models of dark 
energy, we have (? s ~ 1, which explains why such scalar dark energy does not cluster: its own pressure 
resists gravitational collapse. 

It should be stressed that, until recently, no fundamental scalar field had been observed in Nature. 
The discovery of a Higgs particle at the LHC high energy collider has obviously promoted the status of 
fundamental scalar particles. 

In the next section, we will illustrate the dynamics of dark energy on the example of a scalar field 
evolving with time along its potential. This is often referred to as quintessence. Such scalar fields turn 
out to be extremely light: the only dimensionful parameter in the problem being the Hubble constant 
Ho, their mass is HHq ~ ICC 33 eV (see e.g. (233) below). The exchange of such fields leads to a long 


range force: the range is the inverse of the mass (times he), typically ijj l} i.e. the size of the observable 
universe. This force is therefore similar to gravity and can hardly be disentangled from it. Gravitational 
tests such as the test of the equivalence principle thus apply not to the gravitational force alone but to the 
combination of gravity and this new force. This compels this force associated with dark energy to share 
many properties with gravity. One may thus talk of a gravitational type force. 

The second possibility is to modify the left-hand side of Einstein’s equations, i.e. to look for a 


modification of gravity (we already mentionned that possibility in Section 4.1 for explaining observa¬ 
tional data traditionally accounted fo by dark matter). This is a notoriously difficult task because the 
current theory of gravity, general relativity, has passed many stringent experimental and observational 
tests: equivalence principle, Lorentz invariance,.... Any alternative theory should first pass these tests 
equally successfully before being further considered. 

Again, the distinction with the previous case where one introduces a new dynamical component 
is not as clear-cut as it would first seem. Indeed, through field redefinitions, alternate theories of grav¬ 
ity may be rewritten as Einstein gravity plus a dynamical scalar field. Or if one goes to more spatial 
dimensions, the graviton of the higher-dimensional theory may be regarded as a standard 4-dimensional 
graviton plus a collection of scalar (spin zero) or vector (spin one) fields. 

Hence, the distinction between what is often presented as the two ways to account for the observed 
acceleration of the expansion is not so clear. 


5.6 The example of quintessence 

A scalar field d which has reached the minimum do of its potential energy V id) amounts to a cosmo¬ 
logical constant in the form of vacuum energy: its kinetic energy is vanishing and its potential energy is 
the constant V(do)- A more dynamical candidate for dark energy is a scalar field which is still slowly 


evolving in its potential [81-84]. One often refers to this field as a quintessence field. 


To be more explicit, let us consider the general action which describes a real scalar field 
mally coupled with Einstein gravity. 


f , , _ 

m? 1 

j d^xyf^g 



nnni- 


(203) 


Computing the corresponding energy-momentum tensor, we obtain the pressure and energy density (note 
the parallel with Section [23] where we discussed inflation models) 


Pd = 


(204) 
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p<f> 


(205) 


= \i> 2 + v{<t>) 

where, in the latter, we identify the field kinetic energy 4> 2 /2 and the potential energy V (0). The corre¬ 
sponding equation of motion is, if one neglects the spatial curvature (k ~ 0), 

dV 

(f> + 3H<j>= , (206) 


where, besides the standard terms, one recognizes the friction term 3 Hcj) due to expansion. We deduce, 
as expected, 

= -3+ p<f,) . (207) 


We have for the equation of state parameter 

P<j> W -v(<t>) 

uu = — = -> — 1 

P<f> \(j> 2 + V{(t)) 


(208) 


If the kinetic energy is subdominant (<^> 2 /2 <C V(6)), we clearly obtain -1 < < 0. In any case 

-1 < Wej, < + 1 . 

Let us look in more details at the dynamics of such a quintessence scalar field. For this purpose, it 
is useful to identify scaling solutions which we define as solutions where the (f> energy density scales as 
a power of the cosmic scale factor: 


P<f> 


oc a 


-n<l> 


riff) constant 


(209) 


This will allow us to identify two of the main examples of dynamical potential for quintessence: 
- the exponential potential: 


V((/>) = V 0 e 




- the Ratra-Peebles potential [82| 831, 


M A+a 

v (0) = 77v , a > 0 . 


( 210 ) 


( 211 ) 


As we will see, the interest of such solutions is that they correspond to attractors in the cosmological 
evolution of the scalar field. 


Since it follows from (2091 that (),[,/p ( p = —n^H, we deduce from (20V > a relation between 
and the equation of state parameter w 


W'f, = 


x 

3 


( 212 ) 


Flence the scaling solutions that we look for, exist only in epochs of the cosmological evolution where 
the equation of state parameter may be considered as constant (it could still be constant piecewise). 

Since the dark energy (quintessence) density is expected to emerge from the background (radiation 
or matter) energy density, we consider the evolution of a scalar field (f> with constant parameter ?.%, during 
a phase dominated by a background fluid with equation of state parameter 


w 


B 



Following (30), we have a(t) t 2 / n B (n L 


4 for radiation, 3 for non-relativistic matter,...). 


(213) 
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From ( |204| ) and ( |205[ ), we obtain 

4> 2 = ^-p^ v(<f>) = (i -y)p*- (214) 

Hence fi 2 ~ a _n ^ and f ~ t _n <#>/ n s. We thus distinguish two cases: 
i) = n B 

Clearly this implies w 0 = w B > 0 and the quintessence field 0 cannot be interpreted as the dark 
energy component. We have 

(215) 


with X constant. Then 


4> = 4>o + j ln(i/io) 


V (4>) ~ p<f) ~ a n<,> ~ t 


(216) 


Hence, we find a scaling behavior for the exponential potential (210) in a background such that n B = n 0 
(w B = w<p). The solution of the equation of motion (206) then reads 


= - In 


(Vo A 2 


n c 


and the energy density (205 1 


Since H 2 = (p B + pf )/3 ~ [2 /(n B t)] 2 , 


A \ 2 6 — n B 

12 

P(j> — 




P<j> 


X 2 n B t 2 


n L 


P B + P<t> A2 


(217) 


(218) 


(219) 


Hence pp/p B tends to be constant in this scenario. One calls this property “tracking”. This is obviously 
compatible with our initial assumptions only if X 2 > n B . 

What happens if X 2 < n B ? 

It turns out that the scaling solution corresponds to a totally different regime: the scalar field is 
the dominant contribution to the energy density. We do not have to redo the calculation: it is identical 
to the previous one with the only changes w B -0 w<p or n B -0 (for example, H 2 = (p H + pf )/3 ~ 
[2/(rifit)] 2 : the scalar energy density determines the evolution of the Universe). But then (219) reads 
1 ~ n^/X 2 , i.e. 

A 2 


= -1 + y • 

Thus, if X 2 < 2, the scalar field 0 may provide the dark energy component. 


( 220 ) 


To summarize the two regimes that we have obtained for the exponential potential (210»: 

- if A 2 < n B , the scaling solution has = —1 + A 2 /3 and p^Kpcj, + p B ) ~ 1 (<f is the dominant 
species), 

- if A 2 > n B , the scaling solution has w 0 = w B and P(j,/{p B + pf) ~ n ;j /A 2 (the background 
energy density dominates; the scalar field energy density tracks it). 

ii) yh n B 


Then 


—s^+i 
~ t n B 


and we now have 

V((f>) ~ p^ ~ a _n</> ~ t~ 2n 4>l n B ~ f 2 n B ~ n <t . 


( 221 ) 
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Hence, we find a scaling behaviour for the Ratra-Peebles potential (211 1 in a background characterized 
by n B n$ (or iu B ^ w^). We have 


V td 

rich = —— or wa, 

a + 2 


aw B — 2 
a + 2 


The complete solution of the equation of motion (206 1 is 


/ a(a + 2 ) 2 n B M A+a t 2 \ 01+2 
\ 2 [ 6 (a + 2 ) — n B a] ) 


( 222 ) 


(223) 


As we advertised in the beginning, these scaling solutions correspond to attractors in the cos¬ 
mological evolution of the scalar field. We show in Fig. 25 the full phase diagram for the exponential 
potential with A = 2 during matter domination (n B = 3). The coordinates arc x = <j)/(V6Hm 2 p ) and 
y=VV/{VSHm 2 p ). The stable (spiral) attractor is found at x = y = V / 3/8- 



Fig. 25: Phase space diagram for the exponential potential (12101) with A = 2 and n E 


3 (matter domination) f85) 


Exercise 5-2 : We show in this exercise that the solution ( 223 1 of the Ratra-Peebles potential ( | 211\ 
is also an attractor. 


a) Show that a small perturbation 5(f> satisfies the equation 

5f> + H- ,9 [ 6 (a + 2 ) - n B a\ = 0 . 

n B t n B [a + Z)~t.- 


(224) 


b) Look for a solution of the form df ~ f ! and express 7 in terms ofn B and a. 

c) Assume a > 0 and stand ard v alues of n B . Show that the two solutions obtained in b) decay as 

~ the solution (223) is an attractor. 


Hints: b) 


6 — n B \J~ [3a 2 (3n fl - 2)(6 - n B ) - 12a(n| - 16n fl + 12) - 4(?r| - 36n s + 36)] 

2 n B 2 n s (a + 2 ) 

(225) 

c) The reduced discriminant of the second order polynomial in a which is under the square root is 
simply 288/A > 0. The corresponding roots are then negative for 8 — \/b2 ~ 0.79 < n B < 6 and the 
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term is thus negative for a > 0 and n B in thi s range. The square root contributes as an oscillating term 
and the two solutions corresponding to (225) decay as df ~ tr^~ n B^ 2n B. 


The reader may have noticed that the type of field solutions that we have found is very similar 
to those encountered at the onset of inflation. This is not surprising since we are introducing dynamics 
through the slow roll of a field down its potential. Let us look a little closer at the field evolution in the 
case of the Ratra-Peebles potential ( 21 1[ ). 

We have found the attractor scaling solution [82 831 f oc a nB ^ 2+ct \ p 0 oc a -cm s/( 2 +") in the 


case where the background density dominates. Thus p 0 decreases at a slower rate than the background 
density ( p b oc a~ nB ) and tracks it until it becomes of the same order, at a given value a Q . We thus have: 


m.r 


P<j> 

PB 


n B /(2+a) 


2n B /(2+a) 


(226) 

(227) 


Exercise 5-3 : Compute the time t Q at which p 0 ~ p M in terms of M and m p . Check that cf> at t Q 
does not depend on AT. 


2Q+4) 2q 


q+2 


Hints: pm ~ m 2 p /t 2 and p$ ~ M a + 2 t a + 2 give t Q ~ m P 2 M 2 . 


The corresponding value for the equation of state parameter is given by ( |222 i: 

a( 1 + wb ) 


W ( j > = -1 + 


2 + a 


(228) 


Shortly after <j> has reached for a = a Q a value of order m p , it satisfies the usual slow roll conditions 


(using the notations (88 and (891 introduced in the context of inflation) 


1 f m V* \ 2 TTl 2 V" 

€ = ^ f y J = (a/2)(m p /0) 2 < 1 , p = = a (a + \){m p /f) 2 < 1 . (229) 


Therefore (2281 provides a good approximation to the present value of w ( p. Thus, at the end of the 
matter-dominated era, this field may provide the quintessence component that we are looking for. 

Two features are interesting in this respect. One is that this scaling solution is reached for rather 
general initial conditions, i.e. whether p 0 stalls of the same order or much smaller than the background 
energy density [86]. 

The second is the present value of p. Typically, since in this scenario f is of order m p when the 
quintessence component emerges, we must choose the scale AT in such a way that V(m p ) ~ p c . The 
constraint reads: 

AT ~ (H$m 2 p +a ) 1/(4+Q) . (230) 

We may note that this gives for a = 2, M ~ 10 MeV, not such an atypical scale for high energy physics. 


Exercise 5-4 : In the case of slow roll , the equation of motion (206) simply reads '.MI o = —V'{<f). 
a) Under this assumption, show that 

4vr V' 




H 2 


(Pi + Pi) 


(231) 


where the summation is over all components of the Universe. 
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b) Deduce that, in the case where only matter and dark energy are nonnegligible at present time 
to (Qm + ^cj> = 1)> 

(232) 

to 

Hence slow roll requires that ~ 1. 

Hints: a) Use H = -4 nG N Y^iiPi + Pi)- 


jL 

V' 


However appealing, the quintessence idea is difficult to implement in the context of realistic mod¬ 
els j87j|88j. The main problem lies in the fact that the quintessence field must be extremely weakly 
coupled to ordinary matter. This problem can take several forms: 


• The quintessence field must be very light. If we return to our example of the Ratra-Peebles 
potential (2111, V"{m p ) provides an order of magnitude for the mass-squared of the quintessence com¬ 


ponent: 


rritf, 


M 


M \ l+a/2 


m c 


H 0 ~ KT 33 eV. 


(233) 


using (2301. The exchange of such a field leads to a long-range force: the range is typically £jj 0 , the 
size of the presently observable Universe. Since this force has not been observed yet, this means that 
the field </> must be very weakly coupled to matter in order to comply with the constraints imposed on 
gravitational-type forces by the very stringent tests of the equivalence principle. 


• The quintessence field is presently evolving with time. This may generate a time dependence 
of what we call the constants of Nature. Indeed, it turns out that, in modern particle theories, many of 
the constants of nature have a dynamical origin: they are expressed in terms of quantum fields which 
have settled down to their vacuum values in our present Universe. For example, in the models discussed 
above, it is difficult to find a symmetry that would prevent any coupling of the form 


/3— 


rri, 


flld 


(234) 


to the gauge field kinetic term (or any given power of <j )). Since the quintessence behaviour is associated 
with time-dependent values of the field of order m p , this would generate, in the absence of fine tuning, 
corrections of order one to the gauge coupling] But the time dependence of the fine structure constant 
for example is very strongly constrained 


\6l/ol\ < 5 x 10 17 yr 1 . This yields a limit [87]: 


< nr 


dn P Ho 

(<}>) 


(235) 


where {<f>) is the average over the last 2 x Hr years. Let us recall [90] that the non-constancy of constants 


is not compatible with the principle of local position invariance (i.e. independence on the the location 
in time and space where a non-gravitational experiment is performed) which forms part of the Einstein’s 
equivalence principle. This in turn leads to violations of the universality of the weak equivalence princi¬ 
ple m - 


• We have seen that, in the simplest models, the regime of interest is reached when the quintessence 
field cj) value becomes larger than m p . In this instance, as well as in the context of (single field) chaotic 

26 For example, the quantum electrodynamics Lagrangian reads, with our conventions C = —(1/4 d t )F IJ ' v F^. The extra 
term would lead to an effective electromagnetic coupling coupling (l/4e 2 ff) = (l/4e 2 ) + /3(</>/m p ), hence a time-dependent 
fine structure constant a = e^, s /{hc). 
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inflation, there has been discussions whether this remains in the sub-Planck domain: strictly speaking, 
the answer is yes since what characterizes the Planck domain are energy densities of order m A p , whereas 
here p^ remains much smaller because of the specific form of the potential (see e.g. (210) or (2111). It 
remains that, in such a context, one must take into account all non-renormalisable interactions of order 
{fi/m p ) n compatible with the symmetries. 


All the preceding shows that there is extreme fine tuning in the couplings of the quintessence 
field to ordinary matter, unless they are forbidden by some symmetry. This is somewhat reminiscent of 
the tine tuning associated with the cosmological constant. Let us stress however that the quintessence 
solution does not claim to solve the cosmological constant (vacuum energy) problem described above. 
Quintessence may explain the acceleration of the expansion of the Universe but has nothing to say about 
the cancellation of the bulk of the vacuum energy arising from quantum fluctuations. 


5.7 Back to the cosmological constant 

Let us conclude this Section by reviewing some of the attempts to address the problem of the cosmolog¬ 
ical constant. 


5.7.1 Relaxation mechanisms 

In the days where it was believed that vacuum energy was vanishing, one had to look for a mechanism 
to fully cancel the contribution of order m A p . One naturally advocated mechanisms that relaxed the 
cosmological constant to zero through the equations of motions of some dynamical fields. 

For example, in the context of string models, any dimensionful parameter is expressed in terms 
of the fundamental string scale M s and of vacuum expectation values of scalar fields. The physics of 
the cosmological constant and of its relaxation to a vanishing value would then be associated with the 
dynamics of the corresponding scalar fields. 

However, Steven Weinberg [|92j has constrained the possible mechanisms for the relaxation of the 
cosmological constant by proving the following “no-go” theorem: it is not possible to obtain a vanishing 
cosmological constant as a consequence of the equations of motion of a finite number of fields. Wein¬ 
berg's no-go theorem relies on a series of assumptions: Lorentz invariance, finite number of constant 
fields, possibility of globally redefining these fields... All attempts to propose a relaxation mechanism 
have tried to avoid the conclusions of the theorem by relaxing one of these assumptions. 


5 . 7.2 Anthropic considerations 


The anthropic principle approach can be sketched as follows. We consider regions of spacetime with 
different values of to (time of galaxy formation) and t,\. the time when the cosmological constant starts 
to dominate i.e. when the Universe enters a de Sitter phase of exponential expansion. Clearly galaxy 
formation must precede this phase otherwise no observer (similar to us) would be able to witness it. 
Thus tc < t\. On the other hand, regions with t,\ tc have not yet undergone any de Sitter phase of 
re-acceleration and are thus “phase-space suppressed” compared with regions with t\ ~ tc;. Hence the 
regions favoured have t\ ~ tc and thus p\ ~ pm- 


This was quantified by S. Weinberg [15 93), who obtained the following bound: 


PA < 


yPo(l + zg ) 3 , 


(236) 


where po is the present energy density and zq the redshift corresponding to galaxy formation. Using 
Z(_; = 4.5 as originally chosen by Weinberg [ 15j, one finds 0 < ptjP M < 550. More recent observations 
of a galaxy at z = 8.6 |94| or the existence of dwarf galaxies at 2 ~ 10 [95] give a larger anthropic 
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range: 

0 < pa/P M < 4000 . 


(237) 


5.7.3 Emergent gravity 

The alternative approach is to return to the origin of the vacuum energy problem. We stressed in Sec¬ 


tion 5.2 that this problem arises in the context of a quantum treatment of gravity (both h and G N are 
involved). At present we do not have a fully valid theory of quantum gravity. Presumably, it involves 
as well a quantum version of spacetime. It is probable that, just as our notion of continuous and elastic 
matter is only valid in a large distance approximation, our notion of continuous space and time is also 
only valid at large distance. Of course, it remains to be seen by which "objects" one should replace 
continuous space and time, for distances smaller than the Planck length, and what is the corresponding 
theory. In any case, space and time would be emergent notions, and probably also gravity. For what 
concerns us here, it could be that the solution to the vacuum problem should be searched in this deeper 
context. And maybe dark energy is telling us something about this underlying theory. Let us also note 
that, if spacetime is an emergent notion, then its symmetries arc also emergent: one may expect at some 
level violations of Lorentz invariance for example, which lead to violations of Einstein’s equivalence 
principle. 


5.7.4 Holography 

Until now we have considered gravity as a fundamental force which is on the same footing as the other 
three. However, one aspect of gravity is strikingly different from what we encounter with other interac¬ 
tions: it is the phenomenon of gravitational collapse. As we have seen in Section [3T2j if a quantity E of 
(gravitating) energy is localized in a region of spacetime of size R smaller than the Schwarzschild radius 
defined as: 

G E 

R S = 2 (238) 


it undergoes gravitational collapse. This has been used by some (e.g. [96]) to consider that the high- 
energy (ultraviolet) regime of gravity is classical: before reaching Planckian energies, regions of space- 
time undergo gravitational collapse and turn into black holes, which are classical objects. This may have 
some far reaching consequences for the issues we are dealing with here, especially vacuum energy. 

Indeed, let us return to the considerations that led to the estimate p ~ m p /i 3 (see (1831) for the 


vacuum energy density in the context of quantum field theory. Consider a spherical region of radius R 
and energy density given by (1831: p = m A p . Then the total energy reads: 


4i o 47 t . 

e = y r p = Y mp Rrr ^ 


(239) 


But the system will undergo gravitational collapse when R < R$ that is, using (238» i? < ( Rm p ) /3m 


i.e. R > 1 /m p = l p . Hence, for any volume larger than the elementary cell, on cannot concentrate 
a vacuum energy density p = m p , at least in the case (that we consider here) that vacuum energy is 
gravitating: the system is unstable and undergoes gravitational collapse. The maximal energy density for 
a macroscopic region of size R is (E < R/ (2 G N )) 


Pmax = E/(4:irR 3 /3) = 


8 nG AI R 2 


(240) 


Let us extend these considerations to the whole observable Universe of size i? ~ f7 0 1 : we can only store 
a vacuum energy density 


P < 


3ffp 2 
87 tGat 


= Pc ■ 


(241) 
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Taken at face value, this would mean that the vacuum energy density has the value it has because our 
observable Universe is very large. This cannot be true at all times: otherwise, one can easily check that 
the presence of dark energy can be absorbed in a redefinition of Newton’s constant; up to this redefinition, 
the Universe would behave as if there is no dark energy and thus the recent phase of acceleration of the 
expansion would remain unexplained. If pushed to its full consequences, this leads to a new way of 
considering the quantum evolution of the Universe [97]. 


5.8 Concluding remarks 

The most fascinating aspect of the dark energy problem is the number of fundamental questions it con¬ 
nects with: how did the Universe emerge from a quantum state to become so large and so old? why does 
dark energy emerge so late in the evolution of the Universe? does it relate to the nature of space and 
time as we know them? has spacetime emerged from something else? is general relativity the ultimate 
theory of gravity? if not, are its basic principles violated at some scale? what is a quantum state of the 
Universe? what is the status of an observer in such a Universe? are there multiple universes? are there 
more than four dimensions?... 

One of the reasons is that dark energy appears to be connected with vacuum energy, which is the 
most fundamental issue faced by theorists in fundamental physics, an issue that illustrates the difficulties 
encountered at the interface between general relativity, the present theory of gravity, and the quantum 
theory. In some sense, the situation is reminiscent of the one encountered at the end of the XIX th century, 
where one had two very successful theories, Newtonian gravity and electromagnetism (summarized into 
the Maxwell equations). The Michelson-Morley experiment in 1887 was the experimental observation 
that led Einstein and others to reconsider the foundations. Similarly, both general relativity and the 
quantum theory, the latter described at the level of (non-gravitational) fundamental interactions by the 
Standard Model, are extremely successful theories. Moreover, the recent successes of cosmology have 
shown that our picture of the early Universe based on these two pillars is not simply qualitative but is 
supported quantitatively by increasingly precise observations. Is dark energy the signal of a new era? It 
remains to be seen its exact connection with the issue of vacuum energy. But more importantly, it is at 
present a conceptual difficulty, rather than a clear experimental sign of the inconsistency of the overall 
picture. We are still lacking our Michelson-Morley experiment. 

From this perspective, it is reassuring that we have in front of us in the next decade or so a very 
substantial experimental programme, which includes not only increasingly precise and complete obser¬ 
vational data, but also experiments of many types, that might help us identify the road to follow in order 
to reconsider the foundations of physics. 


Appendices 

A Astrophysical constants and scales 

Constants 

The tradition in astrophysics is to use the CGS system. Whereas there are in some specific cases 
useful quantities to be defined (such as the parsec), centimeter and gram seem hardly relevant. We thus 
use here the international system. Note that 1 kg.nW 3 = 1CU 3 g.cnW 3 , 1 J = 10 7 erg, 1 W = 10' 
erg.s -1 . 

Speed of light: c = 299 792 458 m s _1 

Newtonian gravitational constant: G N = 6.6742 x 1 0“ 11 m 3 kg -1 s - 
= G N rrip/(hc ) = 5.906 x 10” 39 

Fine structure constant: a = e 2 /(Aire^hc) = 7.297 X 10~ 3 = 1/137 
Thomson cross section: crx = 87rr 2 /3 = 0.665 barn = 0.665 x 10~ 28 m 2 
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Boltzmann constant: ks = 1.380 x 10 23 J.K 1 = 8.617 x 10 5 eV.K 1 
Planck constant: h = 1.054 x 10~ 34 J.s 

Typical length scales 

Planck length: l p = sj8i:G N h/c 3 = 8.1 x 10 -35 m 

Classical electron radius: r e = e 2 /{Attcq m e c 2 ) = 2.817 x 10 -15 m 

Solar radius: R@ = 6.9598 x 10 s m 

Astronomical unit (au) = Sun-Earth distance = 1.4960 x 10 11 m 
Parsec (au/arc sec): 1 pc = 3.262 light-year = 3.086 x 10 16 m 
Sun-galactic center distance: 10 kpc 
Milky way galaxy disk radius (luminous matter): 15 kpc 
Presently visible universe: 6 Gpc 

Typical mass scales 

Reduced Planck mass: m p = \Jhc/(8irG N ) = 2.14 x 10 18 GeV/c 2 = 3.81 x 10 -9 kg 

Solai - mass: M & = 1.989 x 10 3 ° kg 

Milky Way galaxy mass: 4 to 10 x 10 11 M 0 

Typical luminosities 

Solai' luminosity: L 0 = 3.85 x 10 33 erg/s = 3.85 x 10 26 W 
Typical densities 

Present mean density of the universe: po ~ Pc = 10 -26 kg.in ~ 3 
Interstellai' medium: 10 -22 kg.in ~ 3 
Sun: pq = 1408 kg/m -3 
Neutron stai" 10 18 kg.in 3 

B General relativity 

In the context of general relativity, one defines the Christoffel symbol or affine connection T p which is 
the analogue of the gauge field (it appears in covariant derivatives). It is defined in terms of the metric 
as: 

TV = ^ g pa [<9 m sw + d u g lia - d a g IJP \ , (B.l) 

where g pa is the inverse metric tensor: g pa g aT = Sf 

In the same way that one defines the field strength by differentiating the gauge field, one introduces 
the Riemann curvature tensor: 

R'fiap = daT'fip - dffi p UQ + T p aa T\p - T p p a T\ a . (B.2) 

By contracting indices, one then defines the Ricci tensor and the curvature scalar R 

R IW = R Q liaP , , R = g pv R^ lv . (B.3) 
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One also uses the Christojfel symbols to define the covariant derivatives: 


Vw = V v V fl = d v V, L - Y>\ w V p , 

V= VvV" = d u V» + T>* vp V p . (B.4) 


Exercise B-l : In the case of the Robertson-Walker metric {19), 


a) compute the non-vanishing Christojfel symbols {B.l), 

b) using the fact that the Ricci tensor associated with the 3-dimensional metric 77 is simply 
Ilj j ( 7 ) = 2k"fij, compute the components of the Ricci tensor and the scalar curvature ( B.3), 

c) deduce the components of the Einstein tensor Gdefined in Q: The components of the Einstein 
tensor now read (see Exercise 2-1): 



fa 2 k \ 


Gtt 

= 3 -jt — 2 ) ’ 

(B.5) 


\a z a z ) 


Gij 

= —'fij ( a 2 + 2 ad + k) , 

(B. 6 ) 


Hints: a) Y l ]t = 51a/a, P 7 = ady tJ , T l ]k = T l jk ( 7 ). 

b) R t t = —3d/a, Rij = (2k + da + 2 d 2 ) 7 ij, R = —6 (k + da + a 2 ) /a 2 . 


C Measure of distances 


Measuring cosmological distances allows to study the geometry of spacetime. Depending on the type of 
observation, one may define several distances. 

First consider a photon travelling in an expanding or contracting Friedmann universe. Its equation 
of motion is fixed by the condition ds 2 = 0 (as in Eq. \32\ of Chapter^. One then defines the proper 
distance as 

! r dr ft® rdf 

d(t) = a(t ) / -=== = o(t) / 7 ^ . (C.l) 


/0 \/l — kr 2 


Using 


f to cdt 

It a(t) 


r° cda _ r° cda _ T 
Ja(t) Ja(t) a2 H J 0 


t a (t') ' 
cdz 


H(z) 


we may extract from C.l\ the proper distance at time to-' 
d(t 0 ) =a 0 

Jo 


dr 


sin 1 r k = +1 


/0 sj 1 kr 2 


= a 0 < r 


k = 0 


(C.2) 


sinh 1 r k = — 1 


= ?h ( 




dz 


[^ m (1 + z ) 3 + + z ) 4 + + z ) 2 + ^ a ] ^ 


where £jj 0 = °Ho '• 


If a photon source of luminosity L (energy per unit time) is placed at a distance rfrom the observer, 
then the energy flux (f> (energy per unit time and unit area) received by the obser\>er is given by 


L L 

47 raQr 2 (l + z) 2 47 t d 2 


(C.3) 


The two powers ofl + z account for the photon energy redshift and the time dilatation between emission 
and observation. The quantity d[ = ao'r (1 + z) is called luminosity distance. 
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If the source is at a redshift z of order one or smaller, the effect of spatial curvature is unimportant 
and we can approximate the integral dr/yj 1 — kr 2 in (C.l I by simply r (i.e. the value for k = 0). 
This equation gives 

f to a$cdt 


a 0 r 


a(t ) 


r a ° a 0 cda „ f a ° 
i Q'Q' J a 


da 


(C.4) 


a [1 - qoH 0 {t - i 0 )] 

where we have used the development (47) with tn 0 =^H 0 /c = Hf 1 . Using H${t — tf) ~ (a — ao)/ao <C 


1 and a = ao/(l + z), we obtain for z« 1 


0 I 1 + Qo . 
aor = 1--—z + 


Thus, the luminosity distance reads, for 2«1, 

i + qo 


d L = Ih 0 z 1 - 


z + • • ■ ) (1 + z) — IhqZ ( 1 H- ^—z + ■ ■ ■ 


(C.5) 


(C.6) 


2 J \ 2 

Hence measurement of deviations to the Hubble law (di = £h 0 z ) at moderate redshift allow to measure 


the combination f2 M /2 — 12 a (see (45)). 

Another distance is defined in cases where one measures the angular diameter 5 of a source in 
the sky. If D is the diameter of the source, then D / 5 would be the distance of the source in Euclidean 
geometry. In a universe with a Robertson-Walker metric, it turns out to be a(t)r = oyr /(1 + z). This 
defines the angular diameter distance d ,\ 

dL 


d A = 


(l + zf 


(C.l) 


Several distance measurements tend to point towards an evolution of the present universe domi¬ 
nated by the cosmological constant contribution \ and thus a late acceleration of its expansion, as we 
will now see. 

Exercise C-l : We compute exactly the luminosity distance d[ J = ao'rfl + z) or angular distance 
d:/\ = aor/ (1 + z) in the case of a matter-dominated universe. Defining 

sin -1 r k = +1 

Ck(r) = { r k = 0 , (C.8) 

sinh -1 r k = —1 


use (C.2 1 which reads, in the case of a matter-dominated universe, 

RZ 


to prove Mattig’s formula 


f z dz 

“° a(r)= M K,( 1+2)S+ (i- !i „)(i +2 )t /2 

^M Z + (^M _ 2 ) [a / 1 + ^M Z ~ l] 


a 0 r = 2£ Ho - 




(C.9) 


(C.10) 


Hints: For k f 0, change to the coordinate u 2 = kif l — 1)/ [12(1 + z)\ in order to compute the 
integral (C.9). Using the last of equations (40), which reads ( ,2 Ho / c/,q = fc(12 — 1), one obtains 


C k(r) = 2 Cfc 


A: (12 — 1) 

12 


- Cfc 


' fc (12 — 1 ) 

(1 + z)Cl 


from which (C.10) can be inferred. 


27 at least when analyzed in the framework of the model discussed in this section, i.e. including non-relativistic matter, 
radiation and a cosmological constant. 
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D Perturbations with scalar fields 

We study in this Appendix the perturbations of a scalar field coupled to gravity, following Ref / |99| /. This 
has obvious implications for the study of inflation or dark energy models. 

We consider the most general local action for a scalar field coupled to Einstein gravity: 


m T 


where we have defined 


S = ~ ~Y I d A x^/gR+ / d A Xy/g p{X,(j>) 


1 


X = . 


One may describe this system as a perfect fluid, with the standard energy-momentum tensor (21 1: 

Tfii/ — PPpv W (p p)UjjU u , 

Indeed, varying with respect to the metric, we find 


(D.l) 


(D.2) 


(D.3) 


55 = 


V9 


Xp,xU tl U w - ^pg^u 


5g^ 


= 2 / T i' J d' W 


(D.4) 


with Ufj, = Thus, the energy-momentum tensor has the form (D.3) with the function 

p(X, 0), i.e. the scalar Lagrangian, as the pressure (hence the notation) and the energy density: 


p = 2Xp x - P ■ 


(D.5) 


In the case where p = X — V (0), one recovers ( 205\204 ). 

In the following, a quantity will play an important role; it is the speed of sound: 

c 2 = dP = P, X _ P + P 


dp P,x 2 Xp t x 


(D.6) 


We start with a background metric described by (T79l) (for simplicity, we assume that space is flat: 


k = 0; for the general case, see Ref. [99]) and with a background scalar configuration pit) which 
satisfies (|29p (or equivalently the scalar field equation of motion). 


Perturbing this background, we write in the longitudinal gauge /100 28 
ds 2 = (1 + 2<£>)dt 2 — (1 — 2&)a 2 (t)yijdx l dx 2 , 
where <I> is the Newtonian potential, and we take for the scalar field 

x) = p(t) + 5<p(t, x) . 

Then 


5G° n = 2 


5G°i = 2 


- 3F<t> - 3H 2 ® 
a z 


(D.7) 

(D.8) 

(D.9) 

(D.10) 


s We use the fact that the spatial part of the energy-momentum tensor is diagonal. Otherwise, two different functions $ wand 


would appear respectively as (1 + 2$) in the time component and (1 — 2'P) in the space component /100 J. 
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(DU) 


where, as usual, II 2 = 8nG N p/3. As for the variation of the energy-momentum tensor, we have 

5T° 0 = 5p = p^xSX + p^Sp , 6T°i = (p + p)5Ui . 

Note for the latter that Uq =1, U t = 0 but SUi = (5(p/(p) t i 0. For the former, we use p = 

—3 H{jp + p) = p.xX + p i( pf. One finds 

XrrO _ OTT/ , P + p 

oT o — — 3H{p + p )——|-«— 

S T», = (■> + ,) (^ 

Wfe obtain from dG^ = 8'kG n 8T IJ:V 




(D.12) 

(D.13) 



(a$)' 


( 1+ 4^G„ 0 2(p + /,) A )*’ 
4vrG iV a(p + p) . 


(D.14) 

(D.15) 


The other Einstein’s equations are redundant. We may now define the new variables £ and (V 

o$ = 47rG^, ^ = (D.16) 

(pH a 

which satisfy the equations of motion 

n.(n - 4 - 

(D.17) 


i _ a(p + p) 

? H 2 ( > 


c = 


C 2 77 2 


Defining 


z = 


a 3 (p + p) 

a(p + p) 1 / 2 

^77 


(D.18) 


(D.19) 


and differentiating with respect to conformal time ( f = df/drj = afi and so on), we may write the 
system of differential equations simply as 


1 


? = cpC , C' = , 

z 1 

which can be turned into a single differential equation for (. Indeed, defining v = z(, we find 

z " 

v" — clAv — —v = 0 . 

z 

This can be derived from the following action: 


s = \ 


v' 2 + civ Av + — v 2 


dqdfix . 


(DAO) 


(D.21) 


(D.22) 


If we look for plane wave solutions o.f( |5a71) i.e. v = v^e ' kx , we find two regimes: 

- at long wavelength (k 2 c 2 <C \z"/z\), a non-decaying solution v^ oc z. 

- at short wavelength (k 2 c 2 -C | z"/z\), an oscillating solution oc exp {ikc s rj). 
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Quantization of scalar field in curved spacetime 


We now turn to the quantization of the scalar degrees of freedom. The fact that we are in a non¬ 
trivial background gravitational field brings some new features but, since the background is only time 
dependent, the quantization procedure may be broadly inspired by the flat spacetime case. 

Let us consider a generic scalar field (this could be for example the gravitational potential). The 
standard commutation relations 


[$(r?, x), $( 77 , x')] = [ 11 ( 77 , x), 11 ( 77 , x')] = 0 , [$( 77 , x), 11 ( 77 , x ')] = ^ 3 ( x - x ') (D.23) 


involves the canonical momentum II = 5C/S8^ (we are using here the conformal time 77 ). 

One may decompose the operator <f> over the complete orthonormal basis of the eigenfunctions 
of the Laplace operator. In the spatially flat case that we are considering here, these are simply the 
plane waves: ,\'k( ? ?)e _ * k x (we note that we are making full use of spatial translation invariance, which 
remains a symmetry). We thus write 


$(?7, x ) = - 77 F 


d 3 k 


V2 J (2vr)3 / 2 L 


e- ,k ' x Xk(r7)a[ + e !k ^(r7)a k 


t _i_ c *k.x *, 


where the operators and al satisfy the commutation rules 


[®k> ®k'] — 


a L a l 


= 0, 


®k) 0^7 


= <5 3 (k-k') 


This is consistent with ( D.23) under the condition 

Xk(v)xUv) - XkMXk (v) = 2 * 

The Xk ( r l ) modes satisfy an equation of the type 

Xk(77) + E lxUv) = 0 , 


(D.24) 


(D.25) 


(D.26) 


(D.27) 


where includes a mass-squared term and possibly other contributions such as the one that would 
arise from a non-minimal coupling of the scalar field to gravity. 


In Minkowski spacetime, on constructs a Fock space of states obtained by applying a product of 
creation (negative frequency) operators on the vacuum state | 0 ), defined as the state annihilated by all 
positive frequency operators ak-' fik|0) = 0. This relies on the invariance under the Poincare group 
which gives an absolute meaning to these notions. More precisely, in Minkowski spacetime, the operator 
d/dt is a Killing vector orthogonal to the spacelike hypersurfaces t = constant and the plane wave 
modes e~ lk ' x are eigenfunctions of this Killing vector with eigenvalues —iko = —iu of a given sign. 


In cur\’ed spacetime (see for example the book by Birrell and Davies [1011), the Poincare group 
is no longer a symmetry group of spacetime and correspondingly there is no time-invariant notion of 
positive or negative frequency. There is thus no possibility of agreeing on a specific vacuum state for all 
inertial measuring devices. One may still rely, in some cases, on specific symmetries such as translation 
invariance, conformal symmetry or the de Sitter group to constrain the description of vacuum states. 

[See the review by Mukhanov, Feldman and Brandenberger / |?00 | /] 

We thus choose a given time 770 in order to define a vacuum state |0), ?o such that, for all k, 
a k| 0 )j 7 0 = 0. These annihilation operators are the operator factors of the positive frequency modes 
Xk = X\ i ( r l) ,n th e expansion (D.24 1 (similarly we define the negative frequency modes xfi = Xk(v) )■ 
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If all -Ek are positive, it turns out that one find such modes: they are the solutions of f |D.27[ ) with the 
following initial conditions at time rfy: 


xUm) = E k 1/2 (%), Xkivo) = iEl /2 {m) , 


(D.28) 


consistent with the consistency condition ( D.26\ . Since these solutions obviously depend on r/o, we will 
affect them a superscript (0) in what follows. 

At a later time rj\, we define along the same lines a new vacuum |0) ??1 , which is annihilated by all 
operators b k . These operator s appe ar in an expansion of the type ( D.24\ but with new positive frequency 
modes x k 1 • Since equation (D.27 1 is linear, there is a linear relation between the positive and negative 
frequency modes at r/o and rj\: 

xL 1)+ = a kX'k° )+ + /3 k ,\k° )_ 


Xk } = /3kXk )+ + «kXk } 


|«k| — |/3k I = 1 i 

(D.29) 

where we have used \D.26\ . 

This defines the Bogoliubov coefficients a k and [ff. Obviously we have in parallel for the operators 

K = +a£a k - ^k^k , 

b k = -A k a k + /3k«k • (D.30) 

Let us give an example to illustrate the physical meaning of the Bogoliubov coefficients. We start 
with the vacuum ()), ;o at time r/o and compute at r/i the number of particles b k b k . It is given by 


VO 


(0 


bib k 


9 )r?o — I/3k 


(D.31) 


where we have used ( | D.30 ). Thus, even though we have prepared the system in the vacuum state at 
time r/o, the number of particles is non-vanishing at time rji. Fluctuations can be produced quantum 
mechanically from the vacuum through the coupling of the scalar field to gravity. 

If not all energies are positive, then we cannot define a set of modes through the boundary condi¬ 


tions ( D.28 1. This is in particular the situation encountered in the case of inflation. There, the symmetries 


ofde Sitter space help to define the so-called de Sitter invariant vacuum through the conditions: 
Xk(r?o) = p /2 + ’ Xk0?o) = (3/o +ik- iH' 0 /k) , 


(D.32) 


where FLq = a / (r/o)/a(r/o) = a(r/o)Eo- We recover (D.28) at small wavelength i.e. for k <S TLq. We 
note that, whereas the small wavelength behavior is universal, the large wavelength behaviour strongly 
depends on the choice of vacuum. 


Perturbations 


Let us now apply this formalism to the quantum generation of perturbations. We are interested in 


fluctuations of the Newtonian gravitational potential. Using (D.15) and (D.16), we have 


C = $ 


2 

1 + -- 


P 


2 

+ x 


P 


3 p + p H 


(D.33) 


3 p + p 

Since <I> is constant in any phase where p/p is constant (say matter or radiation domination), then in such 
a phase, ( is simply proportional to the gravitational potential. We thus consider the scalar variable Q 
in what follows and write: 

d 3 k V*k.x Ck(r?)a t+ e ^.x c * (7/)ak 


1 


Ch, x) = ^ 


(27t) 3 / 2 


(D.34) 
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We have Ck = v k/z, where v k satisfies, according to (D. 21), 


V'i + ( C 2 k 2 - 


t’k = 0 . 


(D.35) 


One often characterizes the fluctuations through the power spectrum 'P^(rf) which is defined from the 
2-point correlation function: 

<0|C(x,r ? )C(x + r, 77)|0) = 


'■ fc=+ °° dk sin kr 


' k=0 


k kr 


n(v) 


(D.36) 


Using the decomposition (D.34), one easily obtains 


T>Ct \ _ fc3 \t I 2 - l Uk ! 2 
Vk{l,) 2t r 2 l4kl 2vr 2 |z| 2 


(D.37) 


As we have seen above, Eq. (D.35 1 has two distinct regimes depending of the relative magnitude of c s k 


and z" / z. In the case of slow roll inflation, the main dependence with time in z, as given in ( D.19 1, 
comes from the scale factor. Hence z"/z ~ a"/a ~ ( aH ) 2 . Hence we have to compare c s k with aH i.e. 
the comoving wavelength a/k with the sound horizon length c s /H. 

In the case of short wavelength (smaller than the sound horizon), the normalized solution is 


Vk (2 kcsfl/^ 


ikcsV 


For long wavelengths (larger than the sound horizon), it is 

v k = C k z , 


(D.38) 


(D.39) 


where the constant C k may be obtained by continuity between the two approximate solutions at the scale 
of sound horizon: |Ck| 2 = 1/(2 kc s z s ), with z s the value of z at horizon crossing. 

We thus have the behavior indicated on Fig. [9j the quantum fluctuations are created at small 
wavelength and grow until they cross the (sound) horizon. From then on they grow mechanically until 
they reneter the horizon and are observed in the matter dominated epoch. We thus obtain the power 
spectrum: 


K=^\cJ 


1 H 2 

( H\ 2 

k=aH/c s c sP + P 

w 


k=aH/c s 


The spectral index is defined as 


C 


d lnP k 

ns(k) - 1 = 


din k 

Thus, inflation predicts a departure from a scale invariant spectrum (ns = 1 )■ 


(D.40) 


(D.41) 


Exercise D-l : Show that, in the case of slow roll inflation (Section 2.3), the spectral index is 
simply given by 

ns = 1 — 6e + 2rj . (D.42) 


Hints: In this case, c s = 1. Moreover, 


Use then ( |90| ). 


n s (k) - 1 ~ 


1 dlnP^ 

H dt 


= -6 
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